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Abstract. The progress in the field of topological insulators is impetuous, being 
sustained by a suite of exciting results on three fronts: experiment, theory and 
numerical simulation. Very often, the theoretical characterizations of these materials 
involve advance and abstract techniques from pure mathematics, leading to complex 
predictions which nowadays are tested by direct experimental observations. Many of 
these predictions have been already confirmed. What makes these materials topological 
is the robustness of their key properties against smooth deformations and onset of 
disorder. There is quite an extensive literature discussing the properties of clean 
topological insulators, but the literature on disordered topological insulators is limited. 
This review deals with strongly disordered topological insulators and covers some 
recent applications of a well established analytic theory based on the methods of 
Non-Commutative Geometry (NCG) and developed for the Integer Quantum Hall- 
Effect. Our main goal is to exemplify how this theory can be used to define topological 
invariants in the presence of strong disorder, other than the Chern number, and 
to discuss the physical properties protected by these invariants. Working with two 
explicit 2-dimensional models, one for a Chern insulator and one for a Quantum 
spin-Hall insulator, we first give an in-depth account of the key bulk properties of 
these topological insulators in the clean and disordered regimes. Extensive numerical 
simulations are employed here. A brisk but self-contained presentation of the non- 
commutative theory of the Chern number is given and a novel numerical technique to 
evaluate the non-commutative Chern number is presented. The non-commutative spin- 
Chern number is defined and the analytic theory together with the explicit calculation 
of the topological invariants in the presence of strong disorder are used to explain the 
key bulk properties seen in the numerical experiments presented in the first part of the 
review. 
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1. Introduction 

The field of topological insulators is progressing extremely fast on both theoretical and 
experimental fronts and in the past few years it attracted an unprecedented attention 
from the condensed matter community. This expedited but self-contained review is 
concerned with a less studied aspect of the field, namely the effect of strong disorder 
on topological materials. A material is called topological insulator if it behaves like an 
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insulator when probed deep into the bulk and as a metal when probed near any edge or 
surface cut into the material pQ . This behavior is not triggered by any externally applied 
field and, instead, it is an intrinsic property of the material. For a topological insulator, 
the metallic character of the edges or surfaces is robust against smooth deformations 
of the material, as long as the insulating character is maintained in the bulk of the 
material. Examples of topological insulators are the Chern and Quantum spin-Hall 
(QSH) insulators which will be extensively discussed later. 

The idea behind this review was to bring the attention to a set of analytic 
tools developed for the Integer Quantum Hall-Effect (IQHE) by Bellissard and his 
collaborators in the late 1980's, reviewed in the excellent paper from 1994 [2j. As we 
shall see, these analytic tools can be applied quite directly to many classes of topological 
insulators [3], therefore providing a natural theoretical framework to analytically treat 
the effect of strong disorder in topological materials (for alternative approaches we 
point the reader to Refs. [HE]). Is the work by Bellissard et al relevant to the field of 
topological insulators? We think it is more than ever, at both conceptual and practical 
levels. 

It is important at the conceptual level because the main claim in the field of 
topological insulators is the robustness of the topological properties against disorder. 
This is a huge claim, holding a lot of promises as most of the newly envisioned 
applications are based on it. But to date, we are still missing the hard evidence for 
it, counting experiment, numerics and theory all together. The transport experiments 
on HgTe/CdTe quantum wells j6j [7], the only 2D QSH insulator discovered to date, 
consistently showed a decrease of the conductance with the increase of separation 
between the ohmic contacts. The quantization of the conductance, as predicted by 
theory, was observed only for short contact distances of less than 1 micron (one should 
be aware that HgTe and CdTe materials are extremely difficult to work with so it is hard 
to pinpoint the cause of this behavior). For 3D samples, angle resolved photoemission 
spectroscopy (ARPES) measurements for disordered surfaces [HI [9J HU] seem to indicate 
robust topological extended surface bands, but these conclusions need to be confirmed 
by transport measurements. The transport measurements in 3D topological samples 
have been notoriously difficult [HJ [121 [G2, E] because of a metallic bulk. Recent 
characterizations also showed that band bending near the crystal's surface can trap 
conventional 2D electron states that coexist but also mix with the topological surface 
states [15], making the experiments even harder. But most recently, high quality single 
crystals grown by molecular beam epitaxy have been reported [16j \TT\ UHl [19] (see also 
the important new development in Ref. [20]). When properly doped, these crystals 
display insulating bulk [19] and accurate transport measurements of the surface can 
be recorded. Unfortunately, the first such transport measurements indicate that the 
topological surface states are slightly localized (the localization is believed to be induced 
by inelastic scattering processes). Ref. [19] also reported ARPES measurements, which 
look very similar to the previously reported data [HI |9j [JO] , despite the fact that the 
surface states are localized. The character of the surface states has been also probed by 
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STM measurements [2IJ E2] , which gave evidence that at short scales (the "field of view" 
of the STM is less than 1 /xm 2 ) the states are extended despite the presence of strong 
defects. Shubnikov-de Haas oscillation measurements, which allow one to map the Fermi 
surface (if any), gave an inconclusive picture so far, with few studies [23l EI] reporting 
clear signal coming from a metallic surface and another study [25], done with ultra 
high-quality crystals, reporting no significant contributions from the surface, implying 
an extremely low surface conductivity. The conclusion here is that the experimental 
measurements are converging to a point where the robustness against disorder can be 
rigorously tested but we don't have yet the experimental confirmation of this property. 

The number of existing numerical simulations for disordered topological insulators 
is quite small compared with the number of simulations done in the 80 's and 90 's for 
the Anderson localization. Until recently, there was only one short numerical study [26] 
on the robustness of the bulk extended states in QSH insulators. This study used the 
transfer matrix analysis and was completed for small quasi-one dimensional samples. 
Recently, the transfer matrix analysis was repeated for much larger systems in Ref. [27J, 
re-confirming the existence of robust extended bulk states against disorder in QSH 
insulators. This study seems to contain the most accurate computations to date. The 
transfer matrix analysis was also adopted in [281 [29] for a representative network model 
in the symplectic class, with emphasis on the critical exponents at the mobility edges. 
The properties of the bulk states in a 2D QSH insulator were also probed in Ref. [30] by 
computations of certain topological invariant in the presence of disorder (the systems 
considered in this study are extremely small) . Same method was adopted in Ref. [HI] for 
a 3D QSH insulator with disorder (probably the first simulation in 3D; also very small). 
Ref. [32] presented a level statistics analysis for Chern insulators and computations of 
the Chern number in the presence of disorder. The robustness of the edge states against 
disorder was studied in Refs. [33j I34J |35] by computing the Landauer conductance for 
a disordered ribbon connected to clean leads. All these three studies worked with S z 
conserving models and were limited to small systems (the length of the ribbons was about 
200 lattice sites; this number was 10 8 in the transfer matrix computations of Ref. [27J; 
we also want to mentioned that the theory of the edge states in S z conserving models 
with disorder is firmly established [36]). But despite all these fine numerical simulations, 
we are still lacking systematic studies that combine more than one method and where 
the numerical convergence is rigorously analyzed. From our experience, we can attest 
that, no matter how elaborate these numerical experiments are, there will always be a 
margin of doubt about the localization-delocalization issues, until an analytic proof will 
be available. 

The non-commutative methods have also a great practical value. As pointed out 
in Ref. [4J, the non-commutative formulas of the invariants can lead to extremely 
efficient numerical algorithms, allowing computations of the invariants in the presence 
of disorder for system sizes that are orders of magnitude larger that what was possible 
with more traditional methods [37l[32j[5j. Developing accurate and efficient methods for 
mapping the phase diagram of disordered topological insulators is extremely important 
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for practical applications especially that, as pointed out in Refs. [33l [MJ [27] , the phase 
diagram of a topological insulator can be strongly deformed by the presence of disorder. 

There are additional reasons for writing this review. By its very nature, the 
field of topological insulators can lead to an unprecedented cross fertilization between 
condensed matter physics and various fields in pure mathematics with tremendous 
potential benefits for both fields. We have already seen applications from classic 
Topology [38J E9l SQl Hj, Chern-Simons Theory [42j, Conformal Field Theory, K- 
Theory [431 144] , Random Matrix Theory and nonlinear a-models [451 146] . One hope 
is that we will see many more contributions of this sort from theoretical condensed 
matter physicists and from pure mathematicians. For this reason, we have tried to 
keep the presentation pedagogical and appealing to a wide audience, especially to 
the undergraduate and graduate students looking for good projects, to the theoretical 
condensed matter researchers who like to compute things explicitly, and to specialists 
in Non-Commutative Geometry looking for exciting applications of their field. The 
targeted audience is quite broad and choosing the style of the exposition was not easy. 
Our final choice might not satisfy all readers, but at least we want to let the readers 
know that a great deal of effort was spent on this issue. 

Our discussion is restricted only to the bulk properties of topological insulators in 
two dimensions. Although the current and broadly accepted definition of a topological 
insulator highlights only the robust metallic character of the edge or surface states, every 
known topological insulator seems to display extended bulk states that are robust against 
disorder. This is an extraordinary behavior, especially in two dimensional models. When 
an edge or a surface is cut in a sample of topological insulator, the emerging edge/surface 
states seem to be connected to these bulk states. In fact, the edge/surface states can 
be viewed as these extended bulk states terminating at boundary. For this reason, 
understanding the bulk and the edge/surface properties of the topological insulators is 
equally important. 

We will present several numerical experiments, involving straightforward 
applications of classic techniques such as level spacing statistics. We will follow the 
standard interpretation of the numerical outputs, which will show a clear difference 
between the behavior of a normal and a topological insulator in the presence of disorder. 
This together with a detailed introduction of two models of topological insulators will 
occupy the first part of the review. We also describe here how to define a robust spin- 
Chern invariant for S z non-conserving models (that is, systems that do not decouple 
into independent copies of Chern insulators). Several interesting questions will emerge, 
which will set the direction for the rest of the review. 

The analytic part of the review presents a brisk account of the Non-Commutative 
Theory of the Chern number developed by Bellissard, van Elst and Schulz-Baldes [2]. 
We have reworked certain parts to give the exposition a more "calculatoristic" flavor, so 
that condensed matter physicists who like to compute could easily follow the arguments. 
We have complemented the proofs with discussions and remarks, and tried to keep the 
reader informed at all times about where the calculation or the argument is heading and 
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why do we need to go there. We summarize the arguments before each lengthy proof to 
provide more guidance. We decided to collect the important statements in Theorems, 
Lemmas and Propositions, something to the taste of the mathematicians but that could 
easily irritate other people. One reason for why we chose to do so was to alert the reader 
that these statements have a rigorous proof and that they can be applied with absolute 
confidence. Another reason was that, by doing so, we can state in one place the result 
and the conditions when the result is valid. The last part is especially important for 
our subject because our main goal is precisely to find the most general conditions that 
assures the quantization and invariance of the topological numbers. 

The review includes a presentation of a numerical technique to evaluate the Chern 
[32] and spin-Chern numbers in the presence of disorder. This technique steams directly 
from NCG Theory and allows one to compute the invariants for finite systems without 
imposing the traditional twisted boundary conditions. The finite size formulas converge 
exponentially fast to the thermodynamic limit, given by the NCG Theory. The technique 
allows us to compute the Chern and spin-Chern numbers for large lattice systems and 
large number of disorder configurations (at least one order of magnitude over what is 
currently available in the published literature). 

The review has a section devoted entirely to applications. The class of Chern 
insulators was chosen as the "control case," because they closely resemble the Integer 
Quantum Hall-Effect, already extremely well understood. In this case, the NCG Theory 
gives a full account of all the effects seen in the numerical experiments on Chern 
insulators, presented at the beginning of the review. Calculations of the Chern number 
will allow us to witness explicitly its quantization when the Fermi level is located 
inside the localized part of the energy spectrum, and the failure of such quantization 
when the Fermi level is located inside the delocalized part of the energy spectrum. 
For Quantum spin-Hall insulators, we define the non-commutative spin-Chern number, 
following Ref. [47J, and discuss the conditions when its quantization occurs. Explicit 
calculations of the spin-Chern number indicate again quantization when the Fermi level 
is located in the localized part of the energy spectrum. 

2. Topological insulators: A brief account 

This will be a brief account, indeed. The reason we kept it short is because there are 
now several reviews surveying the evolution of the field and its current status, from both 
theoretical and experimental point of views [HI H9j EH [5]] . Nevertheless, through this 
brief account we want to let the readers know about the impetuous advances that are 
happening right now in the field of topological insulators. 

It is probably a good idea to start from the Integer Quantum Hall Effect, (IQHE) 
which is now extremely well understood. Discovered at the beginning of the 1980's 
[52] . IQHE revealed a truly spectacular manifestation of ordinary matter, displayed in 
the quantization of the Hall conductance and the emergence of dissipationless charge 
currents flowing around the edges of any finite IQHE sample. The intellectual activity 
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spurred by this effect has led to some of the greatest leaps in condensed matter theory. 
Working with a clean periodic system and using Kubo's formula for the Hall conductance 
Off, Thouless, Kohmoto, Nightingale and den Nijs made the famous connection between 
<jh and a topological invariant now known as the TKNN invariant. Using general 
charge-pumping arguments, the Hall conductance was also linked to the classic Chern 
number (see Avron in Physics Today 2003). But it was already clear from the early 
works j53j Ell ESI EE] that impurity states are essential for explaining the Hall plateaus 
seen in the IQHE experiments. The quest for an analytic theory of IQHE that 
includes the disorder has led Bellissard, van Elst and Schulz-Baldes [2] to one of the 
most amazing applications of a new and exciting branch of mathematics called Non- 
Commutative Geometry [57] . This work gives an explicit optimal condition that assures 
the quantization and invariance of the bulk Hall-conductance in the presence of strong 
disorder. Further homotopy arguments for quantization and invariance of the bulk 
Hall-conductance were developed in Refs. [58] and [59j. 

The progress about the edge physics of IQHE satarted with the works by Hatsugai, 
who established in 1993 a fundamental result [60l EI] saying that the number of 
conducting edge channels, forming in an energy gap of the Landau Hamiltonian due 
to the presence of an edge, is equal to the total Chern number of the Landau levels 
below that gap. The technique developed by Hatsugai can deal only with clean systems, 
rational magnetic flux and homogeneous edges with Dirichlet boundary conditions. It 
was only about 10 years later when, using the methods of Non-Commutative Geometry, 
Kellendonk, Richter and Schulz-Baldes established [62], EHl El] a new link between the 
bulk and edge theory, which ultimately allowed them to generalize Hatsugai's statement 
to cases with weak random potentials, irrational magnetic flux and general boundary 
conditions. The equality between bulk and edge Hall conductance was also demonstrated 
by Elbau and Graf [65], soon after the publication of Ref. [63], this time using more 
traditional methods. Further progress was made in Ref. [66], which treated continuous 
magnetic Schrodinger operators and potentials that can assume quite general forms, 
in particular, they can include strong disorder. A similar result was established for 
discrete Schrodinger operators in Ref. [67] . We mentioned that some of these ideas were 
formalized in an abstract setting in Ref. [68j and applications to the edge states problem 
in topological insulators were given in Refs. [691136] . 

The IQHE can be observed only in the presence of an externally applied magnetic 
field. In 1988, Haldane presented a model of a condensed matter phase that exhibits 
IQHE without the need of a macroscopic magnetic field [70] . The systems that behave 
like the one described by Haldane are now called Chern insulators. The time-reversal 
symmetry in these systems is broken like in the IQHE, but it is broken by the presence 
of a net magnetic moment in each unit cell rather than by an external magnetic field, 
as it is the case for IQHE. As we shall see, the techniques developed for IQHE can 
be directly applied to Chern insulators, whose bulk and edge physics [36] is very well 
understood now. The Chern insulators were never found experimentally, even thought 
there is not one known physical reason for this not to happen one day. The Haldane 
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model truly describes the first topological insulator, but because of the lack of the 
experimental evidence, the field of topological insulators took off many years after the 
work of Haldane. 

The interest for a related phenomenon, the spin-Hall-Effect, was picking up in the 
mid 2000's. The effect was predicted decades ago jZHES], and says that a bar made of a 
semiconductor with strong spin-orbit interaction will display spin-polarized edge states 
when an electric charge current is forced through it. The effect was finally observed in 
2004-2005 [731 174] . The search for a quantized version of the spin-Hall-Effect started 
at the time when the results on the classical spin-Hall-Effect were making the news 
[Z3 EH E3 ESI EB EQl ED E21 E3J, but at that time nobody could imagine that there 
are samples displaying a quantized spin-Hall-Effect in the absence of externally applied 
fields. This changed after the discovery of graphene jHH [851 ESI E2], which inspired 
Kane and Mele to propose an explicit model [38] displaying topological edge modes 
carrying a net spin current around the edges. The emergence of the spin-carrying edge 
states is triggered solely by the intrinsic spin-orbit interaction and the flow of the spin 
current is protected by the time reversal symmetry of the model. All materials that 
are not magnetically ordered display the time-reversal symmetry, and there is a large 
number of materials with strong spin-orbit interaction that are not magnetically ordered. 
Therefore, the chance of observing the Quantum spin-Hall-Effect in real materials is 
quite high. The materials exhibiting this effect are now called Quantum spin-Hall (QSH) 
insulators and their hallmark is a dissipation-less spin current flowing along the edges 
of the samples, an effect due to the non-trivial topological properties of the bulk [43] . 

The original calculations suggested that the newly discovered graphene could be 
a QSH insulator. Unfortunately, the spin-orbit interaction is very weak in graphene 
and that makes the experimental detection of the effect very difficult. Nevertheless, 
the race for the discovery of the first QSH insulator was on. HgTe/CdTe quantum 
wells were predicted to display QSH-Effect in 2006 [88j and confirmed experimentally 
in the following year [6j. The first QSH insulator was discovered and a new field 
emerged [1], that of topological insulators defined as materials that are insulators 
in the bulk but metallic along any edge or surface that is cut into the material. 
Unfortunately, the HgTe quantum wells remain the only two dimensional QSH insulators 
discovered to date. In three dimensions, the list of confirmed QSH insulators is quite 
impressive |89l [90l EH E21 [93l EH [95l [96l [971 EHl ESI EH D2H HD21 [Ml EDI EDH [T05] 
and the experimental characterization of these materials is vigorously underway 

una una una una ma urn na na ma na ma cm ma ma usa eb (see also 

the references cited in the Introduction) . Additional classes of topological insulators are 
expected to emerge in the future [1221 11231 1124j . 
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Figure 1. The honeycomb lattice, its unit cell and generating vectors, together with 
some notations used in the main text. 

3. Introduction to Chern Insulators 

3.1. Chern insulators in the clean limit 

A Chern insulator is a periodic band insulator with broken time reversal symmetry, 
with the distinct property of having a net charge current flowing around the edges of 
any finite sample. The time reversal symmetry is broken not by an externally applied 
magnetic field, but by some intrinsic property of the material, such as the occurrence of 
a net magnetic moment in each unit cell. In the following, we will use an explicit model 
to exemplify some of the most important features of these materials. 

The first model of a Chern insulator was introduced by Haldane in 1988 [70], who 
worked with the honeycomb lattice shown in Fig. [TJ The honeycomb lattice can be 
viewed as a triangular lattice with two sites per unit cell (see the shaded region in Fig.[T]). 
The two sites of the unit cell will be labeled by a = ±1 as in Fig. [TJ The triangular 
lattice is generated by the vectors e±\. An additional vector e = — {e\ + e_i) is shown 
in Fig. [TJ which will play a certain role later. The Haldane model involves spinless 
electrons and assumes only one quantum state (orbital) per site, denoted by |n). The 
linear combinations of these states generate a Hilbert space % ) which is equipped with 
the scalar product (n\m) = S nrn . The system is assumed half-filled, which means there 
is one electron per unit cell . 

In its simplest form, the Haldane's Hamiltonian reads: 

#o Cher "= E \n){m\+ £ {( n \n)(m\+e n \rn)(n\}, (1) 

(rim) ((rim)) 

where £n = \{t + irj)a n with a n being the ±1 label attached to each site n, depending 
on how n is positioned in the unit cell (see Fig. [T]), also known as the isospin in the 
condensed matter community. The single (double) angular parenthesis indicate that the 
sum over n runs over all the lattice sites while the sum over m is restricted to the first 
(second) near neighboring sites to n. Notice that a second neighbor hopping always 
connect sites with same a. If we view the honeycomb lattice as a triangular lattice with 
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Figure 2. (a) The bulk spectrum of Haldane Hamiltonian Eq.[T](t = and rj = 0.1) as 
function of (fci, A^)- (b) The energy spectrum of the same Hamiltonian when restricted 
on an infinitely long ribbon with open boundary conditions at the two edges. The 
spectrum is represented as function of k parallel to the ribbon's edges, (c) The local 
density of states (see Eq. |6| of the ribbon, plotted as an intensity map in the plane 
of energy (vertical axis) and unit cell number along the red line shown in panel (d) 
(horizontal axis). Blue/red colors corresponds to low/high values, (d) Illustration of 
the ribbon used in the calculations shown in panels (b, c) and (f, g). The ribbon was 
50 unit cells wide, (e-g) Same as (a-c) but for t = 0.1 and r = 0. 



two sites per unit cell, then the Hamiltonian takes the form: 

(2) 

+ {ia\n,ai)(n + e 1 ,a\ + \n + e 7 , a)(n, a|}, 

where now n denotes the position of the unit cell in the triangular lattice and a = ±1 
is the isospin labeling the two sites of a unit cell. The variable 7 takes the values and 
±1, and e 7 are the vectors shown in Fig. [TJ We actually prefer this later form of the 
Hamiltonian, which will be used from now on. The Hamiltonian depends on the two 
parameters (£,77). We will omit the label "Chern" and use the simplified notation H 
for the Hamiltonian of Eq. [T] throughout the current section. 

In the absence of disorder, we can perform the Bloch decomposition using the 
isometric transformation U from T~L into a continuum direct sum of C 2 spaces: 

U:n^ fc6T C\ U\n, a) = i fe€T e-*»&, (3) 



CONTENTS 



11 



1 


\* 


Normal 
Insulator 


Chern 
Insulator 


C=i 

\K/6 


t 


C=o 









Figure 3. The region of the parameter space where the model of Eq. [I] is in the Chern 
insulating phase and displays the topological edge bands. 



where T is the Brillouin torus T=[0, 2n] x [0, 2n] and: 

6 =fiV «-.=(?Y (4) 





Under this transformation, UH U 1 = ^ ker Hk with: 




i cos(fce 7 ) — rysin(fce 7 ) e ^k-e 1 

H k = y a _. h i ( 5 ) 

e ^ -e 7 — £ cos(fce 7 ) + 7] sin(fce 7 ) 

12 ^ 12 

We denote by e fe ' the two eigenvalues of i/^. The plot of e fe ' as function of fc will be 
referred to as the bulk band structure of the model. 

We now imagine the following numerical experiment. We let the computer pick 
random points in the (t, 77) plane and then perform a computation of the energy spectrum 
for an infinite sample (the bulk spectrum), a computation of the energy spectrum for 
a ribbon shaped sample and a computation of the local density of states (LDOS) for 
the ribbon. The experiment will reveal that, with probability one, the system is an 
insulator because the occupied states are separated by a finite energy gap from the 
un-occupied states, as it is exemplified in panels (a) and (e) of Fig. [2j The bulk band 
spectrum will not reveal major differences between various regions of the (£, 77) parameter 
plane. However, the calculations for the ribbon geometry will bring major qualitative 
differences. For some values such as £=0.1 and 77=0, the energy spectrum for the ribbon 
geometry displays an insulating energy gap, while for values like t—0 and 77=0.1 it 
doesn't. Things become even more intriguing if we look at this spectrum as function 
of the momentum parallel to the direction of the ribbon. Examining panels (b) and (f) 
of Fig. [2| we see that, when t—0 and 77=0.1, the spectrum displays two solitary energy 
bands crossing the bulk insulating gap. For t=0.1 and 77=0, we can still see two solitary 
bands but they don't cross the bulk insulating gap. If we let the computer run for a 
while, picking random points in the (£, 77) plane, it will slowly reveal that this plane 
splits into regions were the model displays bands that cross the insulating gap like in 
Fig. [2](b) and region where the insulating gap remains open like in Fig. [2](f). These 
regions are shown in Fig. [3j 

It is instructive to also take a look at the maps of the local density of states (LDOS): 

p(e,n) = ±Im{(# - e - i0+)" 1 (n, n)}, (6) 
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which will reveal the spatial distribution of the quantum states. The p(e, n) written 
above depends on 3 variables, the energy plus the two spatial coordinates, but for a 
homogeneous ribbon, like the one shown in Fig. [2](d), p(e, n) is independent of the 
coordinate parallel to the edge. Hence p is only a function of energy and one lattice 
coordinate, chosen to be along the red line of Fig. [2](d), in which case we can display 
p using an intensity map. Such maps are shown in Figs. [2](c) and (g). Here, one can 
see that, if the spatial coordinate is away from the edges of the ribbon, there are clear 
regions of practically zero density of states, regions that are perfectly aligned with the 
bulk gaps seen in the band spectra of Figs. [2](a) and (e). When the spatial coordinate 
approaches the edges of the ribbon, the LDOS inside the bulk gap starts to pick up 
appreciable values in Fig. [2](c). This part of the LDOS can come only from the two 
bands crossing the bulk insulating gap in Fig. [5](b). In other words, the quantum states 
associated with these two bands are localized near the edges of the ribbon and, for 
this reason, they are called edge bands. Since the slope de^jdk of a band gives the 
group velocity of an electron wave-packet generated from that band, we can label the 
edge bands as right and left moving. A more detailed analysis of the LDOS will reveal 
that the right/left moving bands are localized on the lower/upper edges of the ribbon, 
respectively (the correspondence will switch if we change the sign of 77). Of course, 
there is a hybridization between two edge bands and a tiny energy gap is opened at 
the apparent band crossing, but this hybridization becomes exponentially small as the 
width of the ribbon is increased. For the ribbon considered in Fig. [2j this hybridization 
can be practically ignored. In fact, if we keep one edge at the origin and send the other 
edge to infinity, that is, we consider a semi-infinite sample, we will observe just one edge 
band crossing the insulating bulk gap. 

If (£, 77) is in the shaded region of Fig. |3j the ribbon is in a metallic state, while if in 
the non-shaded region the ribbon remains in an insulating state. The edge bands seen in 
Fig. [2](b) are called chiral because they connect the valance and the conduction states. 
Due to this feature and provided the bulk insulating gap remains open, those bands 
will not disappear when the Hamiltonian is deformed by either changing the existing 
coupling constants or by turning on additional interaction terms. For this reason, we 
can say that the metallic state of the ribbon is topologically protected. In the trivial 
case, the bands can totally disappear when additional terms are turned on, and what 
typically happens is that the bands sink into the bulk spectrum. When that happens, 
there will be little or no trace of edge spectrum in the LDOS. 

As shown in Fig. [3J if we pick any point (£, rj) at the boundary of the shaded region 
of Fig. [3j we will find that the bulk insulating gap is reduced to zero. We can also see 
some very distinct features emerging, namely, conic points where the bulk bands touch. 
These singular points are called Dirac points and they are actually at the origin of the 
topological properties of the model. When (£=0,77=0) there are two Dirac points, while 
for any other point of the phase boundary there is just one Dirac point. 
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Figure 4. The bulk spectrum when (£, 77) is located at the boundary between the 
topological and normal insulating phases, (a) corresponds to (t = 0,7? = 0); (b) 
corresponds to (t = cos | , 77 = sin |). 

T/ie Chern number 

We give here a brief and formal introduction of the Chern number. Let P denote the 
projector onto the occupied spectrum: 

p = h §d z ~ ^dz, (7) 
where C is a contour in the complex plane surrounding the occupied energy spectrum. 
Under the Bloch transformation, P decomposes in a direct sum of projectors: UPU~ l = 
0fc G 7-Pfc, where is a finite matrix acting on C 2 . It is analytic of fc, except when (£, 77) 
is on the phase boundary. The Chern number is given by the following formula: 

C = ^- J T tT{P k [d kl P k ,d k2 P k }}cPk, (8) 

where a tr" means trace over the two dimensional C 2 space. The integrand in Eq. [8] is 
called the adiabatic curvature and the integral of Eq. [8] can be shown to take only integer 
values, provided the family of projectors are smooth of k over the entire Brillouin 
torus. A plot of the adiabatic curvature is shown in Fig. [5] for the topological phase 
£=0, 77=0.1 (panel a), and for the trivial phase £=0.1, 77=0 (panel b). The plot was 
generated by computing Pj^ on a mesh-grid of 150x150 points and by approximating 
the /^-derivatives by the second-lowest order finite difference. Fig. [5] shows a distinct 
behavior of the curvature when the topological and trivial phases of the Haldane model 
are compared. In both cases, the curvature peaks near the split Dirac points, but in the 
topological phase the peaks have same signature, while in the normal phase the peaks 
have opposite signatures. Consequently, the curvature integrates to a non-zero value 
for the topological case, which is precisely 1, and to for the normal case (plus/minus 
a small numerical error for both cases). A direct calculation will reveal that C takes 
the value inside the trivial phase and the ±1 values inside the topological phases as 
shown in Fig. |3j 

3.3. Chern insulators with disorder 

In this section, we present several numerical experiments on the bulk of Chern insulators. 
They will reveal one of the flagship properties of these materials, manifested in the 
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Figure 5. Plots of the adiabatic curvature as function of fc, when (a) (t,rj) is in the 
topological insulator region and (b) normal insulator region. The top plots show the 
bulk spectrum and they have been aligned with the plots of the curvature in order 
to show that the peaks seen in the curvature plots occur near those k points where 
the bulk bands tend to touch each other. The spikes become less prominent if the 
insulating gap is increased. 

existence of spectral energy regions that contain delocalized states, even in the presence 
of strong disorder. Recall that we are dealing with a 2-dimensional model where, in 
general, the quantum states are localized in the presence of disorder |125j . We will work 
with the following random Hamiltonian: 

H u = H + W J2n,a Vn\n, a) (n, a\ , (9) 

where u) n are randomly distributed amplitudes taking values in the interval [— |, |]. We 
can think of the index uj in H u as the collection of all cj n 's, which in turn can be regarded 
as a point in an infinite dimensional configuration space Q. 

In the following experiments, we used a random number generator to build the H u 
of Eq. |on a lattice containing 30x30 unit cells. Periodic conditions were imposed at 
the boundaries of the lattice. We diagonalized H u and placed its eigenvalues {e l u } i= i^^.. 
on a vertical axis, repeating the calculation 10 3 times, every time updating the random 
potential. The result is a sequence of 10 3 vertical sets containing the eigenvalues {e^} 
for each run, as illustrated in Fig. [6] for different disorder amplitudes W. The level 
statistics was performed in the following way. We picked an arbitrary energy e and, for 
each disorder configuration, we identified the unique and e^ +1 that satisfy: e^<e<e^ +1 . 
Then we computed the level spacings: Ae=e^ +J+1 — e^ +J , letting j take consecutive values 
between — k and k. We have experimented with k— 1-=-5 and the results were virtually 
the same. Fig. [6] was generated with fc=2, in which case, after repeating the procedure 
for all 10 3 disorder configurations, we generated an ensemble of 5 x 10 3 level spacings for 
each e, level spacings that were subsequently normalized by their average. Each diagram 
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Figure 6. (Please rotate at 90 degrees.) Level statistics for the Chern insulator 
(=0.6i (upper panels) and for the normal insulator (=0.6 (lower panels) at disorder 
strengths W=3, 5, 8 and 11. The main panels show the spectrum of for many 
disorder configurations. Level spacings were recorded from a small window around an 
energy e. Shown in light blue is the variance of these level spacings ensembles, when 
e was brushed over the spectrum of H^. The dotted line marks the value of 0.178 
corresponding to the variance of Pque ensemble. The small panels show histograms 
of the level spacings recorded at the marked energies. The histograms are compared 
with the Pgue (blue curve) and Pp i SS on (red curve) distributions. 
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Figure 7. The phase diagram of a Chern insulator as inferred from the numerical 
calculations of Fig. [6j 

to the right of the energy spectrum in Fig. [6] shows the distribution (histogram) of these 
5 x 10 3 level spacings. We picked several values for e and we computed a histogram 
for each value. Imposed over the histograms are two continuous lines, one representing 
the Poisson distribution P(s) = e~ s and another one representing the Wigner surmise 
for Gaussian Unitary Ensemble (GUE), P GUE = ^s 2 e~^ s2 . Imposed over the energy 
spectrum is the variance (s 2 ) -(s) 2 of the 5 x 10 3 level spacings recorded at a large 
(continuous) number of energies. We marked the theoretical value of 0.178 for the 
GUE variance by a dashed line in Fig. |6| The level statistics was performed for several 
disorder amplitudes: W=3, 5, 8, 11, for both the topological phase, t=0, 77=0.6, and for 
the normal phase £=0.6, rj=0. 

Let us focus on the topological case first, shown in panels (a)-(d). Inspecting 
the histograms and the variance in Fig. [6j one can see energy regions where the level 
distribution is Poisson, thus the states in these regions are very likely to be localized. 
But one can also see sharp energy regions where the histograms overlap quite well with 
the 

^Pgue distribution and where the variance converges to the 0.178 value. These regions 
are very likely to contain delocalized states [12611127] . which is quite remarkable since the 
disorder amplitude in all these panels is larger than the bandwidth of the clean energy 
bands seen in Fig. [2} One can also observe in Fig. [6] that, as the disorder amplitude is 
increased, the spectral regions supporting the extended states do not abruptly vanish 
and instead they drift towards each other until they meet and only then they disappear. 
This is the so called levitation and pair annihilation phenomenon, which is a general 
characteristic feature of the extended states carrying a non-zero topological number. 
This will be discussed in depth later. Based on our current observations, the phase 
diagram in the (E F ,W) plane of a Chern insulator with (t,rj) fixed in the topological 
region should look like in Fig. [7| If we examine the normal insulator in FigjHJ panels 
(e)-(h), we see that the spectral regions containing delocalized states are completely 
absent. There is no levitation and annihilation in this case, and instead the extended 
states become localized the moment we turn the disorder on. 

While all that has been said so far about the Chern insulators is just an introduction, 
we already reached the core of our investigation: To establish that the topological 
property of the Chern insulators (to carry a non-zero Chern number) has highly unusual 
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Figure 8. (a) The bulk spectrum of the Kane-Mele model of Eq. [To] (t = 0.1 and 

77 = 0.1 and A = 0.15) as function of (fci,/^)- 3 ) The energy spectrum of the same 
Hamiltonian when restricted to an infinitely long ribbon with open boundary conditions 
at the two edges. The spectrum is represented as function of k parallel to the ribbon's 
edges, (c) The local density of states (see Eq. [6| of the ribbon, plotted as an intensity 
map in the plane of energy (vertical axis) and number of unit cells (horizontal axis) 
along the red line shown in panel (d). Blue/red colors corresponds to low/high values, 
(d) Illustration of the ribbon used in the calculations shown in panels (b-c) and (f-g). 
The ribbon was 100 unit cells wide, (e-g) Same as (a-c) but for (t = 0.4 and 77 = 0.1 
and A = 0.15). 

physical consequences, manifested in the existence of extended bulk states that resist 
localization even in the presence of strong disorder. One of our goals will be to 
demonstrate that the phase diagram in the (E F ,W) plane can be derived analytically, 
using the methods of NCG. 

4. Introduction to Quantum spin-Hall (QSH) insulators 

4-1- QSH insulators in the clean limit 

The first model of a QSH insulator was introduced by Kane and Mele [38], [43] , who 
worked on the same honeycomb lattice of Fig. [TJ but considered also the spin degree of 
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freedom. The Kane-Mele Hamiltonian reads: 



rj-QSH 



E \n,a)(m,a\ 

(rim), a 



+ E a n (t/2 + ir)[S ■ d km x d^ J^n, cr)(m, oj 

((nm)),cr 

+i\ E i e z ■ (S x d nrn )] a ^\n,a)(m,a'\. 

(nm),aa / 



(10) 



Here, S = e x S x + e y S y + e z S z is the spin operator (S x \n,a) = ||n, — a), S^n, a) = 
|a|n, —a) and S^|n,a) = ||n, er)) and a and a 7 the electron spin degrees of freedom, 
taking the values ±1. The Hamiltonian acts on the Hilbert space % spanned by the 
orthonormal basis |n,a), where n is a site of the honeycomb lattice. The simple and 
double angular brackets below the sums in Eq. [TUjhave same meaning as before. Inside 
the second sum, k represents the unique common nearest-neighbor of m and n and d n fc, 
dkm are the displacements shown in Fig.[T] The underline on the vectors in Eq.flOl means 
normalization to unity. The model was built specifically for graphene and the three terms 
in Eq. [TO] represent the nearest neighbor hopping, the second nearest neighbor hopping 
where the intrinsic spin-orbit coupling occurs, and the Rashba potential induced by the 
substrate supporting the graphene sheet or by an externally applied electric field. Half- 
filling will be assumed, that is, two electrons per unit cell, until it is specified otherwise. 
If we view the honeycomb lattice as a triangular lattice with two sites per unit cell, then 
the Hamiltonian takes the form: 



ttQSH _ 



|n, a, cr)(n + crye 7 , —a ) a\ 



+ E {^a,<r\n,a,a){n + e 7 ,«,cr| + + e 7 , a, a)(n, a, a|} r n \ 



-|A ^ acre 



2 ^ 



n, a, cr)(n + crye 7 , —a, — a|, 



where t; a ^ a — \{t + iarj)a and the rest of the notation was already explained. As before, 
we prefer to work with this form of the Hamiltonian, which we will actually do from 
now on. 

In the absence of disorder, we can perform the Bloch decomposition, given by the 
isometry U from the Hilbert space % into a continuum direct sum of C 4 spaces: 



where 
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We have: 



UH Q U- l = ($H' k , with^ 



fceT 




(14) 
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where was given in Eq. |5| and 

_ e 4-^ o )■ (15) 

As before, 7 takes the values 0, and ±1. We denote by e]f the four eigenvalues of the 
Hamiltonian H' k . The bulk band spectrum contains 4 bands, out of which 2 are occupied 
and 2 un-occupied (assuming a half-filled system). 

The parameter space of the model is 3 dimensional (£,77, A). We will let again 
the computer choose random points in this parameter space and then instruct it to 
repeat the numerical experiments already discussed for the Chern insulators. Such 
experiment will reveal that, with probability one, the bulk system is an insulator (see 
panels (a) and (e) of Fig. [8]). Again, by looking only at the bulk band spectrum, we 
will not be able to distinguish any major qualitative difference between different parts 
of the parameter space (t, 77, A) , but the calculation for the ribbon reveals again major 
qualitative differences. For some values such as (t — 0.1, 77 = 0.1, A = 0.15), the energy 
spectrum for the ribbon geometry displays 4 distinct energy bands that cross the bulk 
insulating gap (see Fig. [8](b)). For other values such as (t = 0.4, r\ = 0.1, A = 0.15), the 
spectrum still displays 4 distinct energy bands but they don't cross the bulk insulating 
gap (see Fig. |8](f)). If we let the computer for a while to pick random points in the 
(i, 77, A) space, it will slowly reveal a distinct region were the model displays bands that 
cross the insulating gap like in Fig. [8](b) when restricted to the ribbon, and another 
region where the ribbon has an insulating gap like in Fig. |8](f). 

The plot of the local density of states shown in Fig. |8](f) reveals that the two 
solitary bands marked by (|) are localized on the top edge and the bands marked by (||) 
are localized on the bottom edge of the ribbon. Since the ribbon was 100 units wide, 
there is practically no hybridization between the bands localized at different edges. The 
fact that each edge supports two bands steams from the time reversal symmetry of the 
model and the half-integer value of the spin. The time reversal operation is implemented 
by the anti-unitary operator: 

T = e-™ S vK, ([#o,T]=0), (16) 

where S y is the y component of the spin and K is the complex conjugation. The fact 
that T 2 = — 1 has a distinct consequence in that if ifj is an eigenvector of a time reversal 
symmetric Hamiltonian, than Ti/j is also an eigenvector that is orthogonal to tj) because: 



^,T^) = (T^,T^) = -{^,T^). (17) 

The conclusion is that the spectrum of a time reversal symmetric spin | Hamiltonian 
is always doubly degenerate, a phenomenon known as Krames' degeneracy. For this 
reason, even when considering a semi-infinite sample with one edge, one will necessarily 
observe pairs of right and left moving bands. The Bloch edge Hamiltonian inherits the 
time-reversal symmetry at k=0 and k=n points. At this k points, the spectrum of the 
Bloch edge Hamiltonian is necessarily doubly degenerate, which means the edge band 
crossings occurring at k=0 and k=n cannot be split by any deformation that preserves 
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Figure 9. The left panel shows the spectrum of Hq SH . Its bands below the gap are 
highly entangled. The right panel shows the spectrum of Pa z P. One can see that the 
occupied states have been disentangled into two widely separated bands with non-zero 
Chern numbers. 

the time-reversal symmetry. Edge band crossings occurring at any other k points can 
and are in general split by such deformations. Now, if the number of pairs of chiral edge 
bands is odd, like in the Kane-Mele model, then a simple exercise will show that one 
cannot open a gap in the edge band spectrum by performing all the allowed splittings 
of the edge band crossings. The situation is different if the number of pairs of chiral 
edge bands is even, in which case a gap can be opened, and generically will open under 
deformations that preserve the time-reversal symmetry. This leads to the celebrated Z 2 
topological classification of the time-reversal invariant insulators introduced by Kane 
and Mele [32]. For our simple model, the conclusion is that the ribbon is in a protected 
metallic state. 

Returning to our specific model, if A is set to zero in Eq. [lOj the spin up and spin 
down sectors are left invariant by the Hamiltonian, which is reduced to a direct sum of 
two copies of Haldane Hamiltonian with ( = |(t + irj)a for a=l and ( = \(t — irj)a 
for a—— 1. Concentrating for a moment on the bottom edge of the ribbon in Fig. [8] 
and recalling our discussion from the previous section, one can see that the right/left 
moving edge bands belong to the a=±l sectors. Therefore, the edge bands generate a 
spin flow, because one band carries a a=l spin in one direction and the other carries 
a a=-l spin but in the opposite direction. When A is turned on, the spin sectors are 
no longer invariant under the action of the Hamiltonian and, as a consequence, the 
edge bands will acquire a finite opposite spin component, but still the picture remains 
practically the same. 

Jf.,2. The spin-Chern number for S z non- conserving models 

Time reversal invariant insulators have trivial Chern number. According to Refs. |128l 
129J, these systems are trivial from the general homotopy point of view. However, if one 
insists on preserving the time reversal symmetry, these insulators still display topological 
properties, as we've already seen. There is quite a variety of approaches when it comes 
to the classification of time reversal symmetric insulators j!30l l39l 1431 1471 11311 11321 11331 

[133 [133 [133 [1311 [I3E1 [ffl 

Here we discuss the topological properties of the QSH insulators using the spin-Chern 
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Figure 10. The left panel shows the energy gap of Hq SH for various values of t and A 
and r] fixed at 0.3. The right panel shows the gap of Pa z P. In the left panel, one can 
see a fine white line where the gap of Hq SH is closed. This line delimitates the QSH 
phase. The gap of Pcr z P is positive inside the QSH phase. 



number |138j . While many of the proposed invariants can be, formally, extended to the 
strong disordered case, and this is done usually by employing finite samples and twisted 
boundary conditions, it is not clear at all that the quantization and invariance of these 
extensions survive in the thermodynamic limit when the system becomes gapless due 
to the strong disorder, like in Fig. [6](b) for example. Note that most of these extensions 
tacitly assume that a small spectral gap remain open. In reality, or better said, in 
the numerical simulations with strong disorder, the typical cases are those in which a 
multitude of eigenvalues cross the Fermi level, from below and above when the boundary 
conditions are twisted (hence one cannot avoid Fermi level crossings by moving the Fermi 
level). In this case, the projector onto the occupied states is no longer smooth when 
twisting the boundary conditions, and ad-hoc numerical solutions must be implemented 
to deal with this fact. This is a serious problem, because, even when the states are 
localized near the Fermi level, the number of states crossing the Fermi level can be very 
large and one needs a criterion that specifies when this can create a problem (in some 
sense this is what NCG gives us). Another weak point of these extensions is that their 
quantization can be probed only numerically and the algorithms are not efficient because 
one needs to repeat the diagonalization of the Hamiltonian of a whole sample for a large 
number of boundary conditions. The spin-Chern number is the only invariant to date 
that can be defined in the presence of strong disorder and directly in the thermodynamic 
limit and this is why we focus here exclusively on it. 

Originally [138J, the spin-Chern invariant was defined for a large but finite, squarely 
shaped sample. Twisted boundary conditions were imposed at the boundaries: 

i;(L, y) = e ik * L m y), ^{x, L) = e 2iS ^(x, 0), (18) 

and the projector Pk x ,k y onto the occupied states was calculated for each (k x , k y ) on the 
Brillouin torus. The spin-Chern number was then computed via Eq.[8j The same original 
work has shown, through impressive numerical calculations, that the spin-Chern number 
remains quantized and invariant when disorder is added, even after the insulating gap 
was completely filled with localized spectrum. 

The present discussion follows an alternative definition of the spin-Chern number 
[17], which is more convenient for analytic calculations. The idea is to split the 
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occupied space into two or more sectors with non-trivial Chern numbers. For the 
Kane-Mele model, one can use the spectral properties of Pa z P to achieve just that, 
where a z \n, a) = a|n, a). Indeed, given the Bloch decomposition Eq. 14, one can easily 
compute the Bloch representation of the projector onto the occupied states and then 
form the matrix P^a z P^. In Fig. [9] we chose (t = 0, 77 = 0.6, A = 0.3) and plotted the 
energy spectrum of H' k and the spectrum of Pk& z Pk. Looking at the energy spectrum, 
one can see the two energy bands below the insulating gap being highly entangled. For 
this reason, no topological invariant can be associated with the individual bands. In 
contradistinction, the bands in the spectrum of Pka z P^ are separated by a seizable gap. 
As long as this gap and the gap in the energy spectrum remain open, a Chern number 
can be associated to each individual bands of P^g z P^. If we denote by P^ the projector 
onto the upper/lower eigenvalue of P^a z P^^ we can compute the corresponding Chern 
numbers via 

C± = 2^ [tr{P±[d kl P±,d k2 P±]}d 2 k. (19) 

Since the total Chern number C_ + C+ is zero, C+ — C_ is an even number and we can 
define the spin-Chern number as the integer: 

C S = \{C + -C). (20) 



In Fig. 10 we present the insulating energy gap of Hq sh and the spectral gap of 
P<j z P for r] fixed at 0.3 and t and A varied over a wide range. One can see in panel 
(a) the insulating energy gap closing along a certain line in the (£, A) plane, line that 
delimitates the QSH phase. As one can see, the spectral gap of P<j z P remains open for 
all A and t values inside the QSH phase. The picture remains for any other value of 77, 
showing that bands of Pa z P are always separated by a finite gap and consequently the 
spin-Chern number is well defined. The spin-Chern number takes the values C s — ±1 
(depending on the sign of 77) in the QSH region of the phase diagram, and C s = in 
the trivial region of the phase diagram. 

As a concluding remark, we mention Ref. [151J where analytic calculations of the 
spin-Chern number were carried out for a S z non-conserving model of a QSH insulator. 
These analytic calculations are interesting because they show that the Pfaffian function 
needed in the computation of the Z 2 invariant [43] and the integrand in Eq. 19 are 
closely related. 



4.3. QSH insulators with disorder 

We consider here the Kane-Mele model with diagonal disorder: 

#5 SH = Hq + W £ w w„|n, a, *){n, a,a\. (21) 

We have repeated the numerical experiments presented in the previous section and the 
results are show in Fig. [TT] The computed histograms are now compared with Wigner 
surmise for the symplectic case: P GSE = ^^s 4 e~^ s2 . The variance of this distribution is 
0.104. For the topological case, shown in panels (a-d), the numerical experiment reveals 
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Figure 11. (Please rotate at 90 degrees.) Level statistics for the QSH insulator £=0, 
?7=0.6 and A=0.3 (upper panels) and for the normal insulator £=0.6, 77=0 and A=0.3 
(lower panels) at disorder strengths iy=3, 5, 8 and 11. The main panels show the 
spectrum of for many disorder configurations. Level spacings were recorded from 
a small window around an energy e. Shown in light blue is the variance of these 
level spacings ensembles, when e was brushed over the spectrum of H u . The dotted 
line marks the value of 0.104 corresponding to the variance of Pqse ensemble. The 
small panels show histograms of the level spacings recorded at the marked energies. 
The histograms are compared with the Pqse (green curve) and Pp isson (red curve) 
distributions. 
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again the existence of energy regions where the histograms of the level spacings overlap 
quite well with P GSE and the numerically computed variance takes values extremely close 
to 0.104. These indicate again the existence of delocalized states [126J, which persist 
even for large disorder amplitudes. The levitation and pair annihilation is still visible 
in the upper panels of Fig. [TT] The regions of extended states are absent in the trivial 
case shown in panels (e-h). The levitation and annihilation is also absent in the trivial 
case. 

There is one distinct difference between the Chern and QSH insulators, regarding 
their bulk properties. The spectral regions containing the extended states are reduced 
to a single point for the Chern insulators, while they remain of finite width for QSH 
insulators. Hence, our numerical experiments imply the phase diagram shown in Fig. [T2| 
for the Kane-Mele model when E F and W are varied. 

We would like to end this section with a discussion of the relation between the 
bulk and edge topological properties of QSH insulators. It is now well established that 
the QSH insulators with odd spin-Chern number do present robust edge modes, while 
those with even spin-Chern number do not. This gives the connection between the spin- 
Chern number and the Z 2 classification of the QSH insulators, which is based on the 
edge physics. For example, numerical experiments indicate that the Kane-Mele model 
displays extended edge states over the entire QSH phase drawn in Fig. [T2} But the story 
does not end here. As we will argue in the following, the spin-Chern number protects a 
set of extended bulk states against disorder, like the ones revealed in Fig. [TTJ regardless 
of its parity. We have now gathered enough numerical evidence to announce here with 
confidence that explicit models with C s =2 |152j or with C s —1 but broken time-reversal 
symmetry do posses robust extended bulk states in the presence of strong disorder, even 
though the edge spectrum displays a mobility gap. 



5. The Chern invarint for disordered systems: The Non-Commutative 
Geometry approach of Bellissard, van Elst and Schulz-Baldes 

The definition of the spin-Chern number given in the previous section is based on the 
Chern number. Consequently, the Non-Commutative Theory of the Chern number is 
relevant for both the Chern and Quantum spin-Hall insulators and will be presented in 
depth in this section. 

We start by setting some basic notations. The symbol \\A\\ will denote the operator 
norm: ||A|| := snp ^/{AW\AW) ^ where the supremum is taken over all vectors of norm 
one in the underlying Hilbert space. We will often make reference to the space of 
bounded operators, B{1~L) ) which is the linear space of all A's for which the operator 
norm is finite. When we use the wording "continuous deformation" of an operator we 
mean variations of that operator that are continuous with respect to the operator norm. 
Several additional classes of operators and norms will be introduced later. 

We restrict the discussion to a Hilbert space % spanned by orthonormal vectors of 
the form: |n, a) ) where n G Z 2 is a site of a 2D lattice and a = 1, . . . , K labels the 
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Figure 12. The phase diagram of a QSH insulator as inferred from the numerical 
calculations of Fig. [TT] 

orbitals associated with a particular site. We denote the projector onto these elementary 
states by 7r n>Q! and the projector onto the quantum states at a site n by 7r n : 

K n ,a \n,a)(n,a\ : 7r n ^ a \n,a)(n,a\. (22) 

The theory will be developed for generic orthogonal projectors II that act on this Hilbert 
space (n* = n and II 2 = II). 

5.1. The Fredholm class and the Index of a Fredholm operator 

The Index is one of the main tools of modern topology and will be heavily invoked in 
the following, so we need to begin with the basics of the Index. Given a bounded linear 
operator F, we define its null space as the linear space of its zero modes: 

Null(F) = {0GH|F0 = O}. (23) 

The Index of the operator F is defined as the difference between the number of its zero 
modes and the number of the zero modes of its conjugate: 

Ind{F} = dim Null(F) - dim Null(F*). (24) 

To have a meaning, at least one of the above null spaces must be of finite dimensionality. 
In fact, if we want the Index to be of any use, we must require that both null spaces 
have finite dimensions. We must also rule out the existence of the so called generalized 
wave- functions which obey Fcj) = 0, as it happens when is in the continuum spectrum 
of F. This can be done by restricting ourselves to operators F for which the range FFL 
is a closed space. So we alreay singled out a very special class of operators, called the 
Fredholm class, defined by the following properites |153j : 

Definition. The Fredholm class contains all bounded operators F with the property 
that dim Null(F) < oc, dim Null(F*) < oc, and FFL and F*H are closed spaces. 

The Index, as defined by Eq. [24| takes finite integer values when evaluated on operators 
F from the Fredholm class and it has no meaning when evaluated on operators outside 
of the Fredholm class. For this reason, we always need to make sure that the operators 
belong to the Fredholm class before we evaluate their Index. This is probably the 
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appropriate moment to introduce another class of operators, which contains all operators 
for which the trace operation makes sense. This class is called the trace class and is 
defined as: 

Definition. The trace class S l contains all compact operators A for which J2 i \i{ < oo, 
where {/x^} are the singular values of A ) that is, the eigenvalues of \J AA* . The functional: 

Plk=E^ (25) 

defines a norm on 5 1 , which becomes a Banach spaces. The trace class is a double ideal 
among the bounded operators, that is, BA and AB are both in the trace class if A is 
in the trace class and B is any bounded operator (not necessarily in the trace class). 

The trace operation is finite when evaluated on operators from the trace class, and it 
can be computed as TiA = 1-^10*)? with {<^} being any orthonormal basis in T~L. 

The sum is independent of the chosen basis if A is in the trace class. For A outside the 
trace class, the sum ^(0i|^4|0i) can diverge, be oscillatory or change with the change 
of basis. For this reason, whenever we plan to compute a trace of an operator, we will 
make sure first that the operator is in the trace class, even if this sometimes can be 
more difficult than computing the trace itself. 

We now turn to the question of how to evaluate the Index. If the action of F is 
known explicitly and is simple enough, the Index can be evaluated by using its very 
definition given in Eq. [24} But this is not the case in a large number of situations, in 
which case we must rely on more computationally friendly methods. One such method 
was derived by Fedosov and will be used here. 

Proposition 5.1 (Fedosov's formula). If, for some finite positive integer n, the 
operators (I — FF*) n and (I — F*F) n are in the trace class, then F is in the Fredholm 
class and its Index can be computed as: 

Ind{F} = Tr(I- FF*) n - Tr (I - F*F) n . (26) 

It is a fact that, if (I — FF*) n and (I — F*F) n are in trace class for some n, then they 
are in the trace class for any other integer larger than n. The computation will lead to 
the same Index, independent of which allowed n we choose to work with. In practice, 
however, one tries to work with the smallest possible value for this n. 

When the right part of Eq. [26] is evaluated, it often leads to explicit formulas that 
involve geometric data. In general, a successful and useful Index calculation, which is 
conditioned by the choice of the operator F and the ability to compute the right hand 
side of Fedosov's formula Eq. [26] (or whatever formula we choose to work with), leads 
to statements of the form: 

Ind{F} = /{...} (27) 

where the integral involves objects with explicit geometric meaning, such as the 
curvature in the case of the Chern number. The equality written in Eq. [27] establishes 
that this integral is an integer, something that might be extremely difficult to see by 
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just looking at the integral itself. In modern Topology, the left hand side of Eq. [27] is 
called the analytic Index, while the right hand side is called the geometric or topologic 
Index. The analytic Index provides the quantization (because from its very definition, 
the Index is an integer) and, as we shall see in a moment, also the topological invariance, 
while the geometric Index provides the geometric interpretation and an explicit way to 
compute the actual value of the Index. This is precisely the philosophy that will guide 
us throughout the present paper. The topological invariance follows from the following 
remarkable property of the Index |153l [154J. 

Proposition 5.2. Suppose that the operator F{X) changes continuously with A and for 
all \ 's the operator F(X) stays in the Fredholm class. In this case, Ind{F(X)} is well 
defined and takes the same integer value for all A ; s. 

In other words, the Index of an operator remains unchanged under continuous 
deformations that keep the operator inside the Fredholm class. This principle not only 
gives a very general way to prove the topological invariance of the geometric Index, but 
also allows one to figure out the precise conditions that assures this invariance. For this 
one has to find out how far can F be deformed and still remain inside the Fredholm 
class. This will be exemplified on explicit models, shortly. We end this section by listing 
some additional properties of the Index [154]. 

Proposition 5.3. Let T, S be any Fredholm operators, C any compact operator and U 
any unitary operator. Then: 

i) Ind{TS} = Ind{T} + Ind{S}. 

ii) Ind{T + C] = Ind{T}. 
Hi) IndiUTU- 1 } = Ind{T} 

5.2. The Chern number as an analytic Index: The translational invariant case. 

Assume the existence of a translational invariant projector II acting on the Hilbert space 
with exponentially decaying matrix elements: 

\(n,a\U\m^)\ < ct. e - 7|n - m| , (28) 

with 7 > fixed. Consider the Bloch transformation U from % into a continuum direct 
sum of C K complex spaces: 

U:H^® keT C K , U\n,a) = ±® keT e- ik ^ a , (29) 

where t; a are the single column matrices with entry 1 in row a and zero everywhere 
else. The transformation takes the projector II into a direct sum of projectors on C K : 
UUU -1 = © fcGT IIfc. The projectors 11^ are finite matrices acting on C K space and 
they are analytic of fc, a property that follows from Eq. [28j The following formula for 
the Chern number is familiar to most of the condensed matter theorists and we already 
mention it when we discussed the Chern insulators: 

C = ^ ^tv{tl k [d kl tl k ,d k2 tl k ]}d 2 k. (30) 
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This formula packs geometric and physical content, since the integrand is precisely the 
adiabatic curvature, derived from the adiabatic connection, both having precise and 
deep physical meanings (see Avron in Physics Today 2003). The following result links 
C with an analytic Index. 

Theorem 5.4. Let u : T~L —> % be the following unitary transformation: 

u\n,a) = |g^j|n,a) (n = (rci, rc 2 ), n = n 1 +in 2: n = \ + f). (31) 
Then, if Eq. 28 holds true, the operator (II^II) © H 1 - is in the Fredholm class and: 

Ind{(UuU) © n^} = -L / tr{U k [d kl U k9 d k2 tl k ]}d 2 k. (32) 

Therefore, the Chern number is quantized and invariant to continuous deformations of 
the projector II, as long as Eq. [26] holds true for some finite 7. 

The unitary transformation u describes the effect of an infinitely thin quantum flux 
threaded through the lattice at the spatial point located at n . Recalling our previous 
discussion, it is quite clear that the Chern number appears in Eq. [8] as the geometric 
Index corresponding to the analytic Index of (II'uII) ©IT 1 (IT 1 = 1 — II is the orthogonal 
complement of II). With this relation established, the topological properties of the 
Chern number follow solely from the extremely general properties of the analytic Index. 
If II is the projector onto the occupied states of an insulator, then the existence of 
the energy gap assures the exponential localization of this projector and we can safely 
conclude that the Chern number of the occupied states is quantized and invariant under 
continuous deformations of the insulator that keep the energy gap open. If the energy 
gap closes and opens again, then (IMI) © H 1 - exits and re-enters the Fredholm class, 
in which case the Chern number can assume a different value. This concludes our 
short exercise showcasing the strength of the general strategy that we next apply to 
the disordered case. The above statement will then appear as a particular case of the 
disordered case, which will be treated in full detail. 



5.3. Differential and integral non- commutative calculus for the disordered case 

A disordered configuration will be labeled by the random variable uj ) which takes 
values in a configuration space Q. We assume that the disordered configurations 
cannot be macroscopically discerned, therefore any laboratory measurement involves 
averaging over the random variable uj. The averaging is done with respect to a 
model dependent probability measure dP(uj) on Vt. For example, for the white noise 
considered in the previous sections, dP(uj) = Yi n duJn (this formal expression can be 
given a rigorous meaning). We also consider homogeneous disorder, that is, various 
translations of a sample generate allowed disordered configurations of the original 
sample. Mathematically, this means that each lattice translation by a vector n induces 
a map t n of Q onto itself and, for consistency, these maps must satisfy: 

t n o t rn = t ri + rn . (33) 
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Moreover, the probability measure must be invariant with respect to these maps, that is, 
if Q is a subset of Q and Q f = then dP(uj) = Jg, dP(u). Formally, we can write 
this as dP(uj) = dP(t n uj). We will also require that any subset Q of Q left invariant by 
the maps t n (n G Z 2 ) have zero measure: j^dP{uj) = 0. In other words, we assume 
that the measure dP(uj) is ergodic with respect to the flow t. 

We now introduce the notion of a covariant family of bounded operators. A 
collection of bounded operators is called a covariant family if: 

T n AJT~ x = A tnUJ) for all u G Vt and n G Z 2 . (34) 

Here T n represents the implementation in the Hilbert space % of the lattice translation 
byn: T n \m,a) = \m + n, a). The covariant families form an algebra, since: 

{A, + and {A^B^ (35) 

are also covariant families if {A^^en and {B^uen are covariant families. To make 
sense of infinite sums and infinite products within this algebra, one defines the following 
norm: 

IKAJcJ = sup||4J. (36) 
weft 

The infinite sums and products can be seen as limits, like for example: 

The definition of the bounded covariant families, the algebraic operations and the norm 
defined above leads to a well defined algebra A, which replaces the algebra of bounded 
translationally invariant operators. 

We are now in position to describe the transition from commutative to non- 
commutative calculus. In the translational invariant case, the Bloch decomposition 
allows us to view the operators as matrices, whose entries are functions defined over 
the Brillouin torus. The existence of this classical manifold allows one to define two 
operations over the space of translational invariant operators: a) the differentiation 
with respect to k and b) the integration with respect to k. Let us first describe how 
this operations extend to the disordered case and then comment on the notion of non- 
commutative Brillouin torus. 

The first observation is that, for the Bloch decomposition UAU~ l = © fcG7 -^4(fc), 
we have 

U-'m^rd^AWp = -i[xj, A], j = 1,2. (38) 

The second observation is that: 

^ / Tr{A(k)}d 2 k = Haw,*, ^Tr Area {A}, (39) 

where Tr Area is taken only over the sites included in a finite area. These two observations 
provides the clue on how to proceed. The partial derivatives of a covariant family {i^j^ 
are the covariant families {djA^}^ (j = 1,2) defined by: 

djAu = -i[xj, A u ], for all uj G Q (40) 
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Given that with probability one, 

lim Area ^ 00 ^Tr Area {A} = f dP(u) Tr^A^o}' (41) 
the integration over the Brillouin torus translates into the non-commutative version: 

T({4J) = fdP(u) Tr{7T ^7r }, (42) 
where 7r = |0, a)(0, a| is the projector onto the quantum states at the origin. To 
ease the notation, we will use tr {A} to denote Tr{7r A7ro}. 

The algebra A together with the differentiation and integration defined above can 
be regarded as a non-commutative Brillouin torus [2]. The terminology becomes clear 
if one remembers that the classic manifolds are in one to one correspondence with the 
commutative algebra of smooth functions defined over them, that is, one can reconstruct 
a manifold if he is given this algebra. The Brillouin torus can be recovered from 
the algebra of translationally invariant operators, whose Bloch decomposition leads to 
ordinary functions defined over the classic Brillouin torus. The algebra of translationally 
invariant operators was enlarged here to that of covariant families of bounded operators. 
No classic manifold can be recovered from this algebra, yet we can define a differential 
calculus for it. 



5.4- The Chern number as a analytic Index: The disordered case. 

Given the non-commutative rules of calculus, we can easily spell out the non- 
commutative version of the Chern number: 

C = 2niT{U UJ [d 1 U UJ ,d 2 U UJ ]}. (43) 

The main question now is if the non-commutative Chern number remains quantized and 
invariant to deformations, and under what exact conditions is this happening? 

Theorem 5.5. [Bellissard et al, 1994 fM/] Consider a covariant family of projectors 

{n w }^' if 

r{|[^,nj| 2 }<oc, (44) 

then (n^^II^) © 11^ is, with probability one, in the Fredholm class. Moreover, with 
probability one, its Index is independent of uj and: 

/nd{(n^ry © n^} = 2™ nry^n^, d 2 nj}. (45) 

The above result gives the quantization of the Chern number in the presence of disorder 
but does not provide yet the invariance, which will be discussed later in the paper. We 
want to point out several things. First, one should notice that the condition written in 
Eq. [44] is extremely simple and transparent. It can be reformulated as: 

A 2 = J dP(u) Tr{7r n^ 2 n^ } < oc (46) 

and now one can see that A can be thought as a localization length associated to 11^. 
Second, the condition written in Eq. [44] is optimal, that is, Eq. [44] provides the most 
general condition that assures the quantization of the non-commutative Chern number. 



The remaining of the section is devoted to proving the Theorem 5.5 
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Lemma 5.6. Let II be an orthogonal projection such that (II — uUu*) 271-1 is in the trace 
class for some integer n. Then (IMI) II ^ is in the Fredholm class and 

Ind{(UuU) n^} = -Tr (II - uUu") 2n -\ (47) 

Proof. We will invoke Fedosov's formula. If we take F = (Hull) © 11^, then / — FF* = 
II — 11^11^*11, which can be put into the form 

n - UuUu^u = n(n - uiiu*) 2 . (48) 

Since II commutes with (II — uHu*) 2 , we have 

(n - uuUu*u) n = n(n - uiiu*) 2n . (49) 

Similarly, / - F*F = II - Uu*UuU and: 

(n - uu*uuu) n = n(n - u"Uu) 2n = n^*(n - uiiu*) 2n u. (50) 

Therefore (/ — FF*) n and (/ — F*F) n are in the trace class and we can use Fedosov's 
method to compute the index: 

Ind{(IMI) © n^} = Tr(n - UuUu*U) n - Tr(II - Uu*UuU) n . (51) 

We will use simple algebraic manipulations to bring the above expression to the desired 
form. For this, let us define: 

& =(«°- o) ,6= -i '' andn=i - ; <52) 

A direct computation will show that: 

IMIu*II 
Uu*UuU 

and consequently 




flUUUil (53) 



(n - uuUu*u) n o 

o (n - nti*iMti) n 



(ft - ftl>ftl>ft) n . (54) 



Inspecting Eq. 51, we can see that the Index can be computed as Tr{G , (IT — UUHUH) n } . 
Using the following two identities: 

ft - fte>fte>ft = -[e>,ft] 2 ft and [U, ft] 2 ft = ft[L>, ft] 2 , (55) 

which can be verified through direct computation, we can write: 

Ind{(IMI) © IT 1 } = -Tv{G[U, ft] 2n - 2 ft(ft - UIlU)fl}. (56) 

Because II commutes with [U, II] 2 and G, we can erase one of the two flanking ft's. It 
is convenient to do that in a symmetric way (also note that U 2 = 1): 

Ind{(MI) © IT 1 } 

= -±Tr{G[U, ft] 2n - 2 ((ft - UIlU)tl + ft(ft - UUU))} (57) 
= -lTr{G[f/,ft] 2n - 2 (e>[e>, ft] ft + ft[ft, U]U)}. 
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We can use the identities [77, fl]U = —U\U,U\ and GJJ = —UG to bring one U all the 
way to the front and then use the cyclic property of the trace to move it all the way to 
the back. After these maneuvers we arrive at: 

Ind{(MI) n^} = \Tt{G[U, U] 2n - 2 ([U, ft] ft + IL[U, U])U}. (58) 

Now, if we recognize that [77, ft]ft + ft [77, ft] is nothing else but [{7, ft], we can write the 
compact expression: 

Ind{(IMI) n^} = ^Tr{G[{7, tif^U}. (59) 

This last expression can be cast in the form 

Ind{(MI) n^} = -|Tr{G(ft - UIIU] 21 '- 1 } 

= _i Tr f ( n -^*) 2n_1 ^ 

2 y -(U-^Uu) 271 ' 1 J (go) 

= -|Tr{(n - ullu*) 271 - 1 + u*(U - uTlu*) 271 -^} 
= -Tt^U-uUu^ 271 - 1 }. □ 

The above Lemma gives a very convenient way to compute the Index. We can 
apply this formula to our projectors {11^} by virtue of the following fact. 

Lemma 5.7. The condition T{\[x ) U u; ]\ 2 }<oo assures that (Lt^ — uH^u*) 3 is, with 
probability one, in the trace class. 

We will give the proof of this Lemma at the end, since it is the most technical part of 
the discussion. We now can complete the main Theorem. We will apply the formula of 
Eq. [47] with n = 2. The invariance of the Index with uj is crucial for the computation. 
It can be demonstrated in the following way. Since the translations act ergodically on 
£7, we need to show that the Indexes of (11^11^) 11^ and {I[ trtUJ u\l trtU j) Ht nL j are the 
same for all n G Z 2 . But due to the covariance of {11^}^ and the invariance of the 
Index relative to the unitary transformations, we have: 

lnd{(U tnU uU tnUJ ) U tnU } = Ind{(n w r>T n IIJ IIJ, (61) 

hence we need to show that 

ind{(n w uiL) nj = ind{(n w r>T n n w ) nj. (62) 

The difference between the two operators appearing above is II W (^ — T^uT n ) and 

7/ T*7/T — V ( m+n ° n+m+np \ / fi o\ 

a l n ai n — l^ rjl ^| m+no | \ n +m+n \ J /im ' V U °^ 



But 1^™°! ~~ decays as ^ for large |m|, which means u — T^uT n is a compact 



operator, and therefore the difference between the two operators appearing in Eq [62 
is a compact operator and consequently the two Indexes are indeed the same (cf. the 
general properties of the Index). 



(65) 
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Now, according to Lemma |S~6"1 we have: 

ind{(n wM n w ) e n^} = -Tr {n w - uii w u*} 3 . (64) 

Because the Index is independent of a;, we can average the right part without breaking 
the equality. If we do that, we have: 

n 

= ~ E / ^(w) TV{T> n (IL - uU^fnM 

n 

= ~ E / -L:{7ro(n t _„ w - T>T n n t _^T>*T n ) 3 7r } 

= - E / dP(u) tr {(IL - r*wT ri ri w T*M*T ri ) 3 }, 

n 

where in the last line we made a change of variable t_ n oj — >> uj and used the invariance 
of dP(uj). One should note that we were able to interchange the integral and the sum 
precisely because (11^ — uH^u*) 3 is in the trace class. Let us use the shorthand u n for 
T^uT n and also point out that: 

U n 7T m = 7l m U n = gg^7T m . (66) 

To ease the notation, we will write n for n+n in the following. At this point, we try to 
compute: 

E tr o{(n w - u n iL0 3 } = E tr {(iL - u n n w u*) 

n n,k,m 

x 7r fe (n w - i^n^u* )7r m (n w - w n n w M n } (67) 

n,k,m 

with: 



A(n, k, m) = (1 - $&|) (l - fcgfca, ) i 1 - pSS) • (68) 
We will use the following remarkable identity: 

A(n, fc, m) 2 r Ki{k\ r m2 mife), (69) 

71 

to continue from Eq. [6T| 

. . . = 2ni (hm 2 - mifc 2 )tr {II a ;7r fc II a ;7r m II a; } 

k.m 

= 2m (hm 2 - mifc 2 )tr {n a; [7r fc , IL] [n m , ILJ} (70) 

k,m 

= -2m tro-tIL [i[xi, IIJ, i[x 2 , LLJ]}. 
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Figure 13. Diagrams used in the proof of identity of Eq. [69] The arrow marking an 
angle tells when the angle is positive or negative: a counterclockwise/clockwise arrow 
means positive/negative values. 



One can also show that as long as T{|[as, II a; ]| 2 }<oc, we can interchange the summation 
and integration in the last line of Eq. [65j In this case, we can use the above computation 
to finally conclude: 

ind{(iL^iL)eni} = 27ri J dP{u>) tr {iL H[x l5 nj,-i^ 2 ,nj}. (71) 



It remains to prove the identity written in Eq. |69] and for this we follow Ref. [58j . 
Let us consider A{— n, fc, m) because this quantity has a geometrical interpretation. To 
see that, we place n, k and m in the complex plane at the vertices of our lattice and 
shift n by (|, |) to account for n , as in Fig. 13, We can immediately see that: 

" ~~ e^ 3 , (72) 



-n m—n 



\n\ \k—n\ ' \k—n\ \m—n\ ' \m— n\ \n\ 



where 0^'s are the angles defined in Fig. 13, A direct calculation will show: 



A{— n, fc, m) = — 2i(sin0i + sin0 2 + sin0 3 ). (73) 

We are going to write the summation over n in Eq. [69] in the following way: 

M-n, fc, m) = + 02 + 03) 

- En^i ~ sin 0i) + En(02 - sin0 2 ) - 2J2n(<p3 ~ sin0 3 ). 
Such decomposition is possible because each of the last three sums are absolutely 
convergent, due to the fact that the 0's behave as l/\n\ for large |n|, and consequently 
(j) — sin0 rsj l/|n| 3 . In contradistinction, the first sum cannot be broken in partial sums 
because the resulting sums will not be convergent. Now, since fa— sin fa is antisymmetric 
with respect to the reflection of n relative to the midpoint between the origin and fc, 
and similarly for fa and fa ) each sum in the second row of Eq. [74] is identically zero. 



Furthermore, as exemplified in Fig. [13j we have the simple fact: 

2tt if n is inside the triangle 
+ fa + fa = { 7r if n is on an edge (75) 
if n is outside the triangle 
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At the end, we can conclude that 

— — A(— n, fc, m) = # of grid points inside the triangle, (76) 

n 

with the points on the edges counted only as \. This number is precisely equal to the 
area of the triangle, that is, \\k1m2— (same as — |(fcim 2 — mife)). This concludes 



the proof of Theorem 5.5, □ 



5.5. Macaev spaces and the Dixmier trace 

We are dealing here with the question of weather the operator (11^ — uH^u* ) 3 is in the 
trace class. This question can be answered in a quite elegant way with the help of the 
structures and operation announced in the title of this section. To begin, it is instructive 
to see how (11^ — uYl^u*) 2 fails to be in the trace class, by applying the following simple 
criterion: Given an operator A : A 2 is in the trace class if and only if: 

£ \(m,/3\A\n,a)\ 2 < oo. (77) 

No such simple criterion exists for probing if A 3 is in S 1 . So let's look at the above sum 
when A is 11^ — uH^u*: 

\(m,l3\n u - tdl^*|n,a)| 2 

n,m,a,(3 

2 (78) 
= E 1-kKS) K" + s,m>,a)| 2 . 

n,s,a,!3 

Then, even if we assume an exponential decay of the projector: \(n + s, /^n^n, a)\ < 
ct.e _7 ' s l, which will make the summation over s convergent, the sum still diverges 
because the front term decays only as \n\~ 2 and that is not enough to make the 
summation over n convergent. A second look at what we just said, reveals that the 
sum of Eq. [78] diverges only logarithmically, suggesting that the sum of the singular 
values J^/xf °f (n w — uTluU*) 2 has a weak logarithmic divergence. 

The above discussion revealed the importance of a new class of operators, namely 
those for which ^ fi 2 diverges logarithmically, or more generally, for which ^ /i? 
diverges logarithmically. These classes are now fairly well understood, especially after 
the work of Macaev. 

Definition. For p G [l,oo), the Macaev space A4 P contain all compact operators for 
which: 

1 N 

lim sup uwv $ < 00 • (79) 
We collect the most important properties of these spaces in the following statements. 
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Proposition 5.8. i) A4 P are two-sided ideals in B(H), that is, AB and BA belongs to 
A4 P if A G MP and B is a bounded operator, 
ii) The functional 



defines a norm on M p . With this norm, the Macaev spaces become Banach spaces. 
Hi) If A belongs to M 2 , then A 3 belongs to the trace class. 

The first statement is useful when one needs to establish that a product of operators 
is in MP and the second statement is important when dealing with infinite sums and 
products of operators. The third statement is particularly important for our discussion 
because it tells us that if one can establish that II W — uH^u* is in M 2 , then he can 
automatically conclude that (11^ — uH^u*) 3 is in the trace class. This is in fact the 
strategy that we follow. The classical trace cannot be extended from S 1 to A"! 1 and the 
construction of a new trace over the M l is paramount for our discussion. As we shall 
see, this trace will provide a computationally manageable criterion to test if H^ — uH^u* 
is in M 2 . 

We begin the discussion of this new trace and we follow the notes from Ref . j!55j (see 
also the exposition in [57]). Let us consider first only the positive compact operators, 
i.e. those with the property that > for all cj) G %. If A is such an operator and 

{Hi} denote its singular values, ordered in decreasing order, we introduce the notation 
J1a{N) = ^2 i<N fa- The task is to use /i^(AT)'s and construct a linear functional on 
M 1 , which retains the main properties of the classical trace. To do so, one extrapolates 
P>a(N) to the real numbers by setting J1a{N + A) equal to (1 — X)Pa(N) + \[Ia(N + 1) 
for A G [0, 1], and then defines: 



Let us fix A for a moment and view ta(A) as a functional that assigns positive numbers 
to the operators A G M l . An important result |155j says that 



for any A,BE M 1 . According to this result, the functional appears to become linear 
when A is taken to infinity, which is exactly what is desired, but there is catch. While 
the whole left hand side of Eq. [82] converges to zero as A — >■ oc, due to nice cancelations, 
the individual terms may not have a well defined limit when A — >> oc. While it is true 
that, for A G M l ) ta(A) is bounded for all values of A, ?a(A) can oscillate like cosine or 
sine functions. To solve this shortcoming, one introduces an equivalent relation on the 
linear space C5[3,oo) of continuous bounded functions, where the function A — >> Ta(A) 
belongs to, by saying that / ~ g if the difference f — g vanishes as A — >> oc. If x is any 
strictly positive linear functional on the quotient space C&[3, oc)/ ~, then 




(80) 




(81) 



t a+b (\) - r A (X) - t b (X) = 0(log(log A)/ log A) 



(82) 



Tr Dix {A} = X (ta) 



(83) 
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defines a trace on the positive operators in .A/f 1 , where ta represents the equivalence class 
of the function ta(A) in C&[3, oc)/ ~. This newly introduced trace can be extended by 
linearity over the entire and this gives the Dixmier trace [156J. 

Proposition 5.9. The functional 

M X 3A^ Tr Dlx {A] G C (84) 



is: i) continuous on M 1 when M 1 is equipped with the norm of Eq. 80; 

n) linear: Tr Dlx {A + B} = Tr Dtx {A} + Tr Dtx {B}; 

m) cyclic: Tr Dlx {AB} = Tr Dtx {BA}, for A G M 1 and B G B(H); 

iv) unitarily invariant: Tr Dlx {U AU -1 } = Tr Dlx {A}, for U unitary and A G M 1 . 

v) positive: if A G M 1 is a positive operator, then Tr Dtx {A} > 0. 

We can see that the Dixmier trace has all the characteristic properties of a trace. 
Computing the actual Dixmier trace may not be always easy, especially because of the 
involvement of the linear functional \- However, there are several simple rules which can 
simplify the computation of the Dixmier trace. For example, if lim^^^/x^A^/log 
exists for a positive operator A, then: 

Tr Dix {^}= lim Ji A (N)/logN. (85) 

TV— ^oo 

Since any operator can be written as a linear combination of four positive operators, 
the above rule can be extended by linearity beyond the class of positive operators. 
Another useful and important observation is that Tr Dix {74} = if A is in the trace 
class. Also, for two positive operators A and B such that A — B is positive, we have 
Tr Dix {A} < Tr Dix {£>}. If Tr Dix {£>} can be explicitly computed, maybe by formula given 
in Eq. |85j then this relation gives an explicit upper bound on Tr Dix {yl}, which in many 
instances is as good as computing the trace itself. 



5. 6. Proof of Lemma 5. 7 



Let us lay down the main steps of our proof, which deviates from the original version 
of Ref. [2j. Working with a truncated version 11^ of 11^, to be defined shortly, it is 
first shown that 11^ — vH^u* are in M 2 ) equivalent to say that (11^ — ull^u*) 2 are 
in M 1 . In the next step we establish that Tr Dix {(II^ — ^II^*) 2 } is, with probability 
one, independent of uj. This is an important observation because then we can average 
over uj and J dP(<^)Tr Dix {(II^ — ^II^*) 2 } can be computed explicitly, the result being 
precisely |T({n cj |a3| 2 II a; }), in the limit R oc. Thus, 11^ — ^n^ix* is in M 2 ) and from 



the statement hi) of Proposition 5.8 we can conclude that (11^ — i/II^*) 3 is in the trace 
class. As a remark, notice that the main identity, found in Ref. [2] and detailed below, 
can be written in the equivalent form: 

iv Dix {|}Kn w ]| 2 } = fr{|i[ ; r,n w ]| 2 }, (86) 

which can be seen as the translation to the present context of a central identity in Non- 
Commutative Geometry |157| . Let us now formulate all these into a concise statement. 
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Proposition 5.10. Consider the truncated projector 11^ obtained by setting to zero all 
matrix elements (n + s, /3|n w |n, a) whenever \s\ exceeds R. Then 11^ — uJl^u* are in 
M 2 and, with probability one, the Dixmier trace Tr Dlx { (11^ — ull^u*) 2 } is independent 
ofuu. Moreover, 

Tr Dix {(n* - uH^u*) 2 } = l$dP{uj) ^{IlJa^lU, (87) 

where Xr = ^2\ s \ <r stt s is the truncated position operator. The equality remains valid 
in the limit R — >• oc ; provided T({| [n^, x] | 2 }) < oc. 

Proof. One has: 

- un%u* = T,\ s \<r En 7rn+ a (n w - ^n w w*)7r n = E| S |<i? (88) 

Because of the particular form of the Q s operators, one can easily compute their absolute 
values \Q 8 \ = \JQ* S Q S '- 

\Qs\=T,r 



n-\-s 



\n+s\ 



n 
\n\ 



which are seen to be smaller than 



Er 



n-\-s 



\n+s\ \n\ 



Consequently: 



\\Q 



s\\M 2 < 



Er 



n-\-s n 

\n+s\ \n\ 



7l r 



M 2 



< OO, 



(89) 
(90) 
(91) 



showing that each Q s is in M 2 ] therefore 11^ — idl^*, which is just a finite sum of Q s 
operators, is also in M 2 . 

We now start the computation of Tr Dix {(II^ — uH^u*) 2 }. One has: 

(itf - «n>*) 2 = Y,\ S \<2R E„ W - «n>*) 2 7r n . (92) 

The Dixmier trace of the terms with s ^ is zero. This can be established by considering 
the unitary transformation (here we assume S\ 7^ 0): 

U\n,a) (-1)1^] |n,a), 
which has the following useful properties: 

W7r n+S = ( — 1) L -1 J 7r n+s and 7r n W _1 = (— l)L s i Jyr n . 
From these properties it follows that: 

w (En Wn* - «n>*) 2 7r n ) w 1 
= -En^ +s (n^-tin>*) 2 7r n 

which implies: 

Tr Dix {w (E„ Wn« - «n>*) 2 7r n ) w- 1 } 
= -Tr Dix {EnWn* - «n>*) 2 7r n } , 

and this can happen only if the above traces are identically zero. At this step it was 
established that the Dixmier traces of (11^ — ^II^*) 2 and of ^ n vr n (II^ — ^II^*) 2 7r n 



(93) 
(94) 

(95) 
(96) 
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coincide. One can go one step further and show that they also coincide with the Dixmier 
trace of ^ n a 7r nja (II^ — ull^;u*) 2 7r nia . Indeed, by using the unitary transformation (here 
we assume 7 7^ 0): 

(— 1)M |n, a), 



W\n,a) 



which has the following useful properties: 



+1 



7T, 



and 7r n a W~ 



-1 _ 



(-i)i 



7J7T, 



ra,cn 



(97) 
(98) 



(99) 



one can establish that: 

Tr Dix {W (£ n 7r n , a ^(II* - un*u*) 2 <ir n , a ) W- 1 } 

= "Tr Dix {£ n 7r n>a+7 (n^ - ixn^*) 2 7r n>a } , 

and this is again possible only if both traces are identically zero. 

We are now in the position to show that the Dixmier trace of (11^ — uH^u*) 2 is 
independent of uj. Since the flow t induced by the translations acts ergodically on Q, it 
is enough to show that the Dixmier trace remains unchanged when u is replaced by t r uj, 
with r some arbitrary element of Z 2 . The invariance of the Dixmier trace to unitary 
transformations leads to: 

TW(n£« - «n*x) 2 } - iWM - «n*«*) 2 } 

= Tr Dix {T r (n« - « r n>;) 2 r;} - Tr Dix {(n« - «n>*) 2 } (100) 

The operator appearing inside the last Dixmier traces can be written in the following 
form: 



E El<™,«|n> + S ,/3)| 2 

s<R,f3 n,a 



n-\-r-\-s 
|n+r+s| 



n-\-r 



\n+r\ 



n-\-s 



\n+s\ \n\ 



(101) 



7T r 



The operators inside the first sum are in the trace class, since their singular values are 
less than: 



n+r+s 



n-\-r 



\n-\-r-\-s\ | ti+t* | 



2 




2 




n-\-s n 






|n+s 72 





(102) 



which behave as |n| -3 as \n\ — > oc, hence the singular values are summable. The Dixmier 
trace of trace class operators is zero and this concluded our argument. 
The last steps of the computation go as follows: 

Tr Dix {(Il£-«n>*) 2 } 

i ( 103 ) 

= Tr Dix <^ E ( n > «l(n„ - Mn^M*) 2 |n, a)7r n;Q 



Note that the last Dixmier trace can be computed via Eq. 85 and that the operators 
inside this last Dixmier trace, corresponding to different cj's, commute with each other. 
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Hence, we can average over co and interchange the average with the Dixmier trace, to 
continue 



... = / dP(u) Tr Di J E(^«l(n^-«n>*) 2 |n,«)7r niQ 

In, a: 

= Tr Dix \^2jdP(u){n,a\(U^ - uU*u*) 2 \n, a)n n A 

(104) 

= Tr Dix \ ^JdP(u)(0,a\T*(U^ - uY[*u*) 2 T n \0,a)^ a 



= Tr Dix {^JdP(uj)(0,a\(Ul nUJ - UrJ^_ nU u* n ) 2 \0,a)ir n>a \ 

At this point we use the invariance of the measure dP{uj) to the flow t to continue as: 
• • • = Tr Dix (e / dP{u)(0, a\(U2 - u n ll^u* n ) 2 \0, a)n n>a 



Tr Dix ^ E / dP(u) E (0,^(11^11^0, a) ^-^^a) (105) 

\s\<R 



= fdP(u) E (o,«|(n w 7r s n w |o,«)Tr Dix (E 2 

\s\<R,a In J 

The remaining Dixmier trace can be explicitly evaluated and the result is | |s| 2 . At this 
point we arrived at: 

Tr Dix {(n^ - un*u*) 2 } = i$dP{uj) E (caKiLkl^iLlo,^, a o6) 

\s\<R V ) 

which is precisely the statement of the Proposition. □ 
The limit R — >> oc can be safely taken if J dP(uj) t r { 11^ | | 2 11^ } < oc, in which 
case the result shows that 11^ — uli^u" is in M 2 class and consequently (11^ — uU^u*) 3 
is in the trace class. 

5. 7. Invariance of the Chern number 

So far, the theory established the optimal conditions that assures the quantization of 
the non-commutative Chern number corresponding to a generic covariant family of 
projectors. We call these conditions optimal because they cannot be improved any 
further. The next question is under what conditions does the Chern number remain 
constant under various deformations of the projectors? The invariance of the Chern 
number with the shift of the Fermi level was already discussed in in the original work 
by Bellissard et al [2j, and was reconsidered in Ref. [58]. The invariance of the Chern 
number under more general deformations was considered in Ref. [59] . These works based 
their analysis on estimates done directly on the Chern number formula. Here we will 
base our analysis on the special properties of the Index. There are two complications. 
Even for smooth deformations of the Hamiltonian, the projector onto the spectrum 
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below Ep does not change smoothly, because eigenvalues will inherently cross the Fermi 
level. Secondly, some of the operators 11^11^ may exit and re-enter the Fredholm class 
during the deformation, because a finite localization length A does not imply that all 
11^11^ are in the Fredholm class (only almost all). 

Theorem 5.11. Consider a deformation {^(r)}^^ of a covariant family of projectors, 
such that: 

(i) T({|[^IUt)]| 2 }) < A 2 < oo, for all r G [0,f]. 

(ii) There exists a covariant family of self- adjoint operators {Q^t)}^^ that are, with 
probability one, in the Fredholm class, and such that 

\\QM-QUr')\\ and \\(QM(r) - (^n w )(r')|| < ct.\r-r'\, (107) 

with the constant independent of uj or r . 

In these conditions, the Chern number o/{n^(r)} wG n is well defined and assumes the 
same integer value for all r G [0, f]. 

Proof Condition i) assures that [II a; (r)^n ct; (r)] 11^ (t)- 1 is in the Fredholm class and 
that Ind{[n w (r)^II w (r)] © II^t)- 1 } is the same for all uj e fi, except a possible set of 
measure zero. Given the general properties of the Fredholm operators, same can be said 
for Q^(r)([n w (r)^n w (r)] ©II w (r) ± )(5 w ( r r) and, since the Index of a self-adjoint operator 
is zero, we have: 

108 

= Ind{Q w (r)([IL(r)wIL(r)] © TL^Q^r)}. 

From the general theory of the Index [153J, we know that each set O n , containing all 
Fredholm operators with Ind{F} = n, is an open set in the space of bounded linear 
operators. The open sets O n are disjoint and one cannot continuously move from one 
O n to another without exiting the Fredholm class. In the following, we will use the 
wording "almost all" to indicate that a statement holds true for all cj's, except for a 
set of measure zero. Now, because of the invariance of the Chern number with uo ) a 
set O n will either contain almost all operators Q^^t^U^^uII^^t)] © Ii UJ {r)^Q UJ {r) ) or 
almost none of them. There is no alternative. This means that if the value of the Chern 
number changes, almost all operators will have to simultaneously migrate from one O n 
to another. This is prohibited because all these operators change continuously with 
r. □ 
The complication with the operators Q u {r) is unavoidable because, in the concrete 
examples, the projectors H u (r) will not change continuously with r, since eigenvalues 
will unavoidably cross the Fermi level during the deformation. Although this section 
dealt with a general family of projectors {II w } wG ^, we will like to switch now to the 
case when 11^ is the spectral projector onto the spectrum below E F \ II^^P^ = 
x(Etu;—E F ), with x{— x ) being the Heaviside function and {H^}^^ a covariant family 
of Hamiltonians. In this case, the operators can be constructed via a standard 
procedure. 
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Corollary 5.12. Let {H^^ea be a covariant family such that 

{n,a\H UJ \m,(3)=0if\n-m\>R, (109) 

with R arbitrarily large but fixed. Consider a deformation {H^t)}^^ such that the 
above relation holds true for r in some interval [0, f ] and such that: 

\{n,a\Hu{r) - Hu{T')\m,f3)\ <ct.\r-r'\, (110) 

for all r and r f G [0,f]. The constant above is assumed independent of ' uj, r or r f . If 

T({[x, x (H u (t) - E F )} 2 }) < A 2 < oo (111) 

for all r G [0, f] and the spectrum of H^r) near E F is regular, that is, with probability 
one, it has finite degeneracy, then the Chern number of {x{H u {r)—E F )} u ^ is well 
defined and assumes the same integer value for all r G [0,f]. 

Before we start the proof, we have the following remark. The condition written in 
Eq. Ill implies spectral localization [2], that is, the spectrum near E F is pure point. In 
the following, when we use the wording "localized spectrum," we refer to that part of 
the spectrum where Eq. Ill holds true, together with the regularity condition of being 
finitely degenerate. 

Proof. We take to be (j){H UJ —E F ) ) where cj){x) is a smooth function equal to one 
everywhere except for a small interval around the origin. We require that (j)(x)>0 
except at x=0, where we require 0(O)=O and d N cj) / dx N \ x=0 =0 for all TV's. Such (j) can 
be constructed via standard procedures. Let (j)±{x) be the smooth functions that are 
identically zero for negative/positive x and equal to cj>{x) for positive/negative x. If E F 
is in the pure point spectrum, then c/)(H U j—Ef)'H is a closed subspace of T~L and since 
the degeneracy of the localized spectrum is finite, dim Null (^(i?^— E F ))<oo. In other 
words, (^{H^—Ep) is in the Fredholm class. Moreover: 

Qu{[Pu,uP u \@P£)Q„ f 

112 

= - E F )v4-{H„ - E F )} </>+(H u - E F ). 

As opposed to the operators (j)±(H cu —E F ) are changing smoothly under the 
deformation and this concludes the argument. 

In the following we demonstrate that cj)±(H UJ —E F ) changes continuously during the 
deformation. To see this, we employ a functional calculus introduced by Dynkin |158j 
(see also Helffer and Sjostrand |_159j ) . We consider a domain V in the complex plane 
defined by all those z—u + iv with \v\ < vq and u in an interval taken large enough so 
that it contains all the spectra of H^r). For any smooth function (p(x), such as (j)(x) 
and (j)±(x), and positive integer N : we can construct a function f:V^C such that: 

(i) /jv(^ 3 z) =( f( z ) when z is on the real axis, 
(h) \dzf N (z,z)\ < a N \lm z\ N . 
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Such / is called an almost analytic extension of ip and one has: 

<p(H u (r)) = ±f v d s f N (z,z)(z-H u (T))- 1 <Pz. (113) 
The matrix elements {n ) a\^{H u {r))—Lp{H UJ {r f ))\rri ) (5) can be computed as: 

± f v d 2 z d 2 f N (z, z) Yl ( n > a \( z ~ H^r))- 1 ^, a') 

P,r,a',(3' (114) 

x (p, a'\5H\p + r, /3'){p + r, f3'\(z - H u {r'))- 1 \m, /5), 

with 5H=H LU (r)—H UJ ( K T f ) and the summation over r restricted to \r\<R. Since the norm 
of an arbitrary bounded operator B satisfies the relation: 

||5|| 2 < supj] \ {n,a\B\m,fi)\ sup ^ |(n, a\B\m, /3)|, ^15) 

m,/3 n,a n,a rn,(3 



we need to show that sums of the matrix elements of Eq. 114 are convergent and 



proportional to |r— r'\. We will use an estimate that follows from a standard application 



of the Thomas-Combes technique |160j . and is valid for all operators satisfying Eq. 109 
[36J. This estimate is: 

\(n,a\(u + iv - H^r))' 1 ^ /3)\ < ct.H^e"*-^ (H6) 

with e small but finite. All constants appearing above are independent of uj and r. We 
then have: 

\{n,a\<p(H u (T))-<p(H u (T'))\m,P)\<ct.\T-T'\ 

X f„ M du f V0 dv \ v \N-2J2 e -e\v\(\n-p\+\ P+ r-m\) 

J Em J-VO 1 1 

P,r ( U7 ) 
< ct \ T _ r /| £0 dy v N-l e -ev{\n-m\+l)/2 

<ct.\r-T , \{\n-m\ + lf- N . 
Taking strictly larger than 5 will insure that the matrix elements are summable. □ 

6. Numerical Evaluation of the Non-Commutative Chern number 

Before we turn to applications, we discuss a practical formula to evaluate the non- 
commutative Chern number in the presence of disorder. The formula was announced in 
Ref. 1 321. The non-commutative Chern number formula: 



C = 2m f dP(a;)tr {P w [ - i[x u P u ], -i[x 2) P u }}} (118) 

makes sense only for infinite samples. Here we will look for a finite size formula that 
converges exponentially fast to the above exact formula. We start from the observation 
that, even in the clean limit where the Bloch decomposition is available, we still cannot 
compute C exactly, and various numerical approximations must be adopted. Now, 
whatever formula we find for the disorder case, it must also work for the clean case. 
Therefore, our first goal will be to find a real-space, finite size formula for C in the clean 
case. We start from Eq. 1301 
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It is clear that the integrand of Eq. [30] can be evaluated only on a discrete set of 



points of the Brillouin torus, most often a finite N x N mesh Ai: k n = (niA,n 2 A), 
77,1,2 = 1, • • • , N, where A = ^. We will use the notation A 1 = (A, 0) and A 2 =(0, A). 
The integral over k is then replaced by a Riemann sum: 

/{...Kfc^£ fcneA ,{...}A 2 . (119) 

Assuming for a moment that the integrand can be computed exactly, such action will 



is 



introduce errors that vanish exponentially with AT, because the integrand in Eq. [30 
analytic of k. 

We now turn our attention to the integrand. While the projector can be indeed 
computed exactly on the mesh considered above, its derivatives with respect to ki^ can 
be computed only approximately, using finite differences: 

dkiPk ~> $kiPk n E!=l Cm[Pk n +mA l A^-mA,], (120) 

with Q < y. We need to find the appropriate combination of the c m coefficients 
so that the finite difference approximation introduces only exponentially small errors. 
Being a smooth function defined on a torus, the projector i\ can be decomposed as an 
absolutely converging sum of the form ^2 neZ 2 b n e lkn . Therefore we need to find a good 
finite difference approximation for functions of the form e lkx , with x arbitrary but fixed. 
We have: 

(d ki -5 ki )e ikx = i(x t - 2Z®= 1 CrnSmmAx i y kx . (121) 

If we choose c=(ci, c 2 , . . . , cq) to be solution of the the following linear set of equations: 

Ac T = |(f,0,0,...,0) T , Av=j*-\ i,je{l,...,Q}, (122) 

then (dk i —Sk i )e lkx =o( K A 2 Q~ 1 ), an exponentially small error if we take Q of the order y. 
Our next task is to translate the discrete fc-space formula: 

C = h E kn tT{P kn [S kl P kn ,8 k2 P kn }}A\ (123) 

into a dual real space formula. Note that Pk n is nothing else but a discrete Bloch 
decomposition of the projector corresponding to the original real space model restricted 
to Sl NxN lattice with periodic boundary conditions. Let Hn be the Hilbert space 
generated by |n, a) with n restricted to this finite lattice. The discrete Bloch 
transformation is given by: 

U : Un -+ ® kneM C K , U\x, a) = ® knEM ±e~ ik ^e a (124) 

and 

U- 1 (® kneM ^(*»)) = E jf E e ik - x r(km)\n,a). (125) 

n,a km^M 

Now note that Eq. |123| can be equivalently written as: 

C = -2niTr{® kneM (P kn [5 kl P kn ,h 2 Pk n })}, (126) 
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Figure 14. Deformations along any of the rays (1), (2) or (3) inherently lead to a 
change of the Chern number from 1 to 0. This allows one to infer the phase diagram 
pictured above. 



where the trace is taken over the whole Hilbert space (&k n eM ^ K • Since the trace is 
invariant under isometries, we can map this trace to the real space: 

C = 2mTr{U-\ 0, eM P kn [S kl P kn , 8 k2 P kn ])U}. (127) 



is: 



The real space representation of the finite difference of Eq. |120 



(128) 



= (x,a\J2Z=i c ™( e ~ im&AiPeim±Ai - e imiAi Pe- im * Ai )|t/,/3). 
which leads us to the desired real-space formula for the Chern number: 

C = E n Jn, a\P{ - i\x x , P\ , -i[x 2 , P\]\n, a), (129) 

where 

[x„ PJ = i YH=i c m {e- im ^Pe im ^ - e ™ A ^ e -^). (130) 



For the clean limit, the formula of Eq. 129 is completely equivalent with the discrete 



fc-space formula Eq. 123 In this case, it is known explicitly and is numerically confirmed 



that this Eq. |129| converges exponentially fast to the exact value of the Chern number. 



For example, computed on a 30x30 (40x40) lattice for the Chern insulator model with 
77=0.6, the formula gives (7=0999998 (0.999999998). In the presence of disorder, we can 
evaluate the Chern number via: 

Cu = -$^ n Jn,a\P„[-ilx u P„\,-ilx 2 ,P UJ \}\n,a). (131) 

When the Fermi level is in a mobility gap, we expect this formula to converge 



exponentially fast to the exact expression given in Eq. |118, as TV is taken to infinity, 



a fact that we indeed observe in our numerical simulations. Note that Eq. |131| is self- 
averaging, so the dependence on uj disappears in the thermodynamic limit. 
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7. Applications 

7.1. Chern Insulators 

The Non-Commutative Theory can be applied to arbitrarily complex Chern insulator 
models. We will, however, return to the simplest Chern insulator model discussed at the 
beginning. As an application, we will derive the phase diagram shown in Fig. [7] using 
only the invariant properties of the Chern number. 

We consider a smooth deformations of the model by changing one parameter that 
can be 1) a constant added to the Hamiltonian, which will effectively change the Fermi 
level Ep, 2) the disorder strength W, or 3) the hopping amplitude ( a . We denote this 
parameter by the generic symbol r, like we did before in the general discussion of Section 



5.7| H u is of course of finite range and the bound: 

\(n,a\H u (r)-H u (r f )\rn^)\ <c*.|t-t'|, (132) 

is uniformly satisfied for r in an arbitrarily large but finite interval [0, f]. In other words, 
we are in the conditions of the Corollary |5.12[ Given all the above, we can assert with 



absolute confidence that the Chern number of the occupied states is an integer as long 
as the Fermi level resides in a localized part of the energy spectrum. The Chern number 
can change its value during the deformation only if the Fermi level crosses a region of 
delocalized spectrum. 

We now focus to what happens when we move the Fermi level or increase the 
disorder strength W. Since all the other parameters are kept fixed, we can identify the 
system with a point in the (E F) W) plane. Let us start from a half-filled system and 
a Hq that is in the Chern insulator part of the clean phase diagram shown in Fig. [3j 
We now consider deformations of the systems along different rays that start from this 



initial point, as shown in Fig. [14} Note that the energy spectrum of H u , represented 
by the gray region in Fig. [l4j becomes wider with W but stays finite. Now, moving 
along any such ray, the Chern number will eventually become zero because one of the 
following reasons: the projectors P^ become the identity (ray 1), all the bulk states 
localize when W crosses a threshold value (ray 2) |161j . or the projectors P^ become 
null (ray 3). Given these facts, a phase diagram like the one discussed in Fig. [7] emerges, 
which necessarily contains a region where (7=1, a surrounding region where (7=0, and 
a region of delocalized states separating the (7=0 and (7=1 regions. The numerics and 
scaling arguments hint that the region of delocalized states is infinitely thin. 

Fig. [15] shows the numerical values of the Chern number, computed via the 
technique presented in Section [6] for seven Fermi levels. The model parameters were 
given the same values as in Fig. [6](a) and the calculations were performed for 30x30 and 
40x40 lattices. In absolute agreement with the general theory and with the discussion 
above, the calculations show that, for Fermi levels that are in between the delocalized 
spectral region, identified here with the region where the covariance is very closed to 
the value of 0.178, the Chern number has the quantized value of one or tends to this 
quantized value as the lattice size is increased. The Chern number decreases to zero 
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-4-202 



Figure 15. The lower panel shows the energy spectrum and variance of the level 
spacings, which already appeared in Fig. [6ja). The upper panels show the numerical 
values of for many disorder configurations, corresponding to the Fermi levels 
indicated by the vertical marks. The blue/red data in these panels correspond to 
computations completed on 30x30/40x40 lattices, respectively. 



as the Fermi level passes the region of delocalized spectrum and after that it assumes 
again a quantized value of zero. 

7.2. Quantum spin-Hall Insulators 

7.2.1. The non- commutative spin-Chern number. We restrict the discussion to the 
model of Eq. [Toj In the clean limit, we have argued that the QSH insulator can be 
characterized by the spin-Chern number C s defined in Eq. [20} The definition involved 
the spectral projectors P ± of the operator Pd z P onto its positive/negative spectrum, 
more precisely the Chern number of these projector. Given that and using the non- 
commutative theory of the Chern number, we are now in the position to extend the 
definition of the spin-Chern number to the disordered case. Let P^ be the projectors 
onto the states below the Fermi level and let P^ be the spectral projectors of P^&zP^ 
onto its positive/negative spectrum, respectively. Let C± be the corresponding Chern 
numbers, defined through the no-commutative formula: 

C ± = 2ni J dP{u>) ti {P£[-i[x u P±],-i[x 2 ,I*]\}. (133) 

Then the non-commutative spin-Chern number is defined as 

C S = \(C + -C). (134) 
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The next Sections will discuss the conditions in which C s is well defined and takes 
quantized values. 

7.2.2. Weak disorder regime. We call weak disorder the regime when the amplitude 
W is small enough such that a gap of size 2Ai is still present in energy spectrum of 
and a gap of size 2A 2 is present in the spectrum of P^dzP^. The Fermi level is assumed 
to be inside the energy gap. The second condition was actually observed to hold true 
in all our numerical experiments, with weak and strong disorder. 

One should not mistakenly think of this weak disorder regime as the regime in 



which a perturbative approach is justifiable. For example, W=?> in Fig. [TT] belongs to 
our weak regime, but note that W=?> is larger than the width of the clean bands. Such 
disorder is quite large and definitely cannot be treated perturbatively. 

Proposition 7.1. In this weak disorder regime, the matrix elements (n, a, a|P^|m, /?, a') 
decay exponentially with the separation \n — m\. Consequently, the spin-Chern number 
C s is well defined and takes a quantized value. Moreover, C s cannot change its value 
during continuous deformations H u (t) that leaves the two gaps open at all times. 

Proof. First of all, let us establish analytically that such weak regime exists. We have 
seen that, in the clean limit, the gap of Po z P remains open for all H Q inside the QSH 



part of the phase diagram (recall Fig. 10). Assume now that H is fixed in the QSH 
part of the phase diagram and let us increase W. As long as the spectral gap of 
remains open, the projectors P^ are continuous of W . At W=0, we know as a fact that 
\\(Pa z P — 0) _1 || w 1, so we can be sure that || (P u d- z P u — 0) _1 || < oc, for W smaller than 
a critical value. This critical value depends on where exactly is H located in the phase 
diagram. The conclusion here is that the origin is outside the spectrum of P^o^P^, at 
least for small W ', and that the spectrum is divided in positive and negative islands and 
we can safely define the projectors P^. 

Our next task is to establish the exponential localization for Pj 1 , which can be 
accomplish via a straightforward application of the Combes-Thomas technique [160J. 
We decided to present a fairly complete discussion for the readers that are not familiar 
with this technique. We do not seek here optimal estimates, which could be obtained 
with the technique developed in Ref. |162j . but rather the simplest proof possible. We 
will first establish the exponential localization of the full projector P u . Consider the 
invertible map: 

|n, a, a) — » U q \n, a, a) = e _<7 ' n |n, a, a), q G R 2 . (135) 

Due to the finite range of P^, one will easily find that H u (q) = U q JiJJ~ x has the form 
(q — \q\) Huj(q) = Hu + qV{q) ) with V{q) independent of uj and uniformly bounded by 
some V if q is restricted to a finite interval [0, q}. For a z located a distance Ai away 
from the spectrum of H Ul we have: 

||2 - H u (q)\\ = \\z -11^ + qV(q)\\ > \\z - HJ - qV = A, - qV. (136) 
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Figure 16. Each box shows the energy spectrum of (left) and next to it the 
spectrum of P^a z P^ (right) for many disorder configurations at different values of W 
and strengths of the Rashba potential. The other parameters were fixed at 77=0.6 and 
t=0. The Fermi level was fixed at the origin in these calculations. 



In other words, z — H UJ (q) is invertible for small q : and 

||(*-ff w ( q ))-i||<^-^. (137) 

Noticing that 

(z - HUq))- 1 = U q (z - Hu)-^, (138) 
which can be verified by applying z — H u; (q) to both sides, we can conclude that: 

ei< n -^\(n,a,a\(z - H^lmJ,*'^ < 1 (139) 

Ai - qV 

Taking q parallel to n—m, we obtain the exponential decay of the resolvent with the 
separation \n — m\. The exponential decay of the projector follows from the formula 



/ dz 

(n,tt,a|P w |m,/5,cr') = j> — {n,a,a\(z - H^)' 1 ]™, (3, a') 



(140) 



where C is a contour surrounding the spectrum below the spectral gap. All points on 
this contour can be taken a Ai distance away from the spectrum of H UJ) so the above 
estimates can be applied without modifications. 
We now turn our attention to P^. Since 

f dz 

(n,a,a\P±\m,/3,cr') = J — {n,a,<r\(z - P w ^P w ) _1 |wi, /3,a'), (141) 
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we can try to repeat the arguments we used for P w , this time working with P w o z P w 
instead of H w . This will work because 

P^OzPJJ- 1 = P u (q)* z P u (q), (142) 

and, since P u (q) is a smooth function of qr, 

PM&zPM = Pu&zPu + qR{q), (143) 

with R(q) uniformly bounded by some R if q is restricted to a finite interval [0, q]. From 
here, we can repeat the arguments used for P w , to conclude: 

c + e -q\n-m\ 

\(n,a,v\P£\m,P,o')\ < — g-. (144) 

A 2 - git 

The exponential localization of P^ can be established for any r of a deformation that 
keeps the two gaps open. Moreover, for such deformations, the projectors P^(r) change 
continuously with r. These two facts prevent any sudden jumps of the spin-Chern 
numbers. □ 

7.2.3. Beyond the weak disorder regime. In all our numerical experiments, we have 
observed a clear separation of the spectrum of P^a z P u into positive and negative parts. 
In fact, we consistently noticed a gap in the spectrum of P u a z P u near the origin, whose 



origins are not clear to us at this moment. We exemplify this in Fig. [16j where we show 
the energy spectrum of and of P^c^P^, when P^ is computed with E F fixed at 0. For 
all values of A and W shown in this figure, we can clearly observe the spectral gap for 
P UJ a z P UJ) even when the energy gap of was completely filled with localized spectrum. 
More importantly, the spectrum of P a a z P u closer to the origin was always seen to be 



localized. In Fig. [TTJ we show the size dependence of the spectrum of P^a z P^ for W=5 
and different Rashba potential strengths. A clear spectral gap near the origin remains 
visible in all the plots and the spectrum closer to the origin was found to be localized. 
Owing to the NCG theory of the Chern number, we can state the following fact. 

Proposition 7.2. The spin-Chern number is well defined and takes quantized values 
as long as the spectrum of P^&zP^ near the origin is localized. The invariance of C s 
can be also established if some additional assumptions are considered, similar to those 
of Refs. [58] and JM/. 

The main issue now is to establish under what conditions is the spectrum of 
Pu&zPuj localized near the origin. The following simple argument will show that the 
localization lengths of the projectors P^ diverge whenever the localization length of the 
full projector P w diverges. This shows that the phase defined by C s —1 is contained 



within the boundaries defined by the metallic phase established in Fig. 11, Here is the 
argument: 

\(n,a,<T\P u \0,a,a)\ 2 = |(n, a, a\P~ + P+|0, a, a)\ 2 

(145) 

< 2(|(n,a,a|Pj|0,a,a)| 2 + |(n,a,a|P+|0,a,a)| 2 ), 
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30 x 30 35 x 35 40 x 40 




Figure 17. The spectrum of P^&zP^ at W=h and various strengths of the Rashba 
potential, computed on lattices of increasing size. The other parameters were fixed at 
?7=0.6, t=0 and the Fermi level was pinned at the origin. 



which gives, 

|A(P w ) 2 <A(P-) 2 + A(P+) 2 . (146) 

Therefore, whenever A(P LJ ) 2 diverges to +00 the localization lengths of P^ (which are 
equal due to the time-reversal symmetry) will be forced to diverge to +00. 

The bulk theory of our QSH insulators will be completed if we can establish that the 
phase defined by C s —1 actually extends all the way to the metallic phase inferred from 
Fig. 11. For this we need to establish that the localization length of the full projector 
P^ diverges whenever the localization lengths of P^- diverge. This will not only show 
that the C s —1 phase extends all the way to the metallic phase, but will also imply the 
existence of such metallic phase surrounding the QSH phase. In other words, we will be 
able to establish the phase diagram of the QSH model in a manner similar to what was 
presented for Chern insulators. This, unfortunately is a highly technical issue and falls 



beyond the scope of this paper. Nevertheless, the NCG statement from Proposition 7.2 



the inequality of Eq. |146 and the numerical technique for evaluating the spin-Chern 
number, have given us the ground to study the remaining issue via numerical means. 

In Fig. 18, we show what happens when we lower the Fermi level and cross the 
spectral region containing the extended states. The bottom panel of the figure shows 
the energy spectrum of at together with the variance of the level spacings. 

The mid panels show the spectrum of P u a z P u with P^ computed for the marked Fermi 
levels. The upper panel shows the numerical value of the spin-Chern number. The 
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Figure 18. The lower panel shows the energy spectrum and variance of the level 
spacings, which already appeared in Fig. 11 'a) (tilted horizontally here). The mid 



panels show the spectrum of P^cr^P^ with the Fermi level pinned at the vertical marks, 
and the variance of the level spacings. The upper panels show the numerical values of 
the spin-Chern number for many disorder configurations, corresponding to the same 
Fermi levels indicated by the vertical marks. The computations was completed on a 
40x40 lattice and for 1000 disorder configurations. 



calculation was performed on a 40x40 lattice, for the QSH model with i=0, 77=0.6 and 
A=0.3. We can see that, when the Fermi level is away from the spectral region containing 
the delocalized states, the spectrum of P^&zP^ displays a gap near the origin, and the 
variance of the level spacings for this spectrum indicates that the spectrum near the 
origin is localized. Similarly, the spin-Chern number takes quantized values whenever 
the Fermi level is away from the energy region of delocalized states. This behavior was 
seen in all our numerical experiments, suggesting that indeed the C s =l phase extends 
all the way to the metallic phase. It is interesting to note that the variance in the 
mid panels collapses entirely onto the value 0.178 when the Fermi level is in the region 
of extending states, suggesting that P^&zP^ becomes a stereotypical Gaussian random 
matrix when P w becomes delocalized. 

8. Conclusions 

The discussion we just concluded dealt with two simple but fundamental questions 
for the field of topological insulators: How to define quantized topological invariants 
in the presence of strong disorder and what are the physical properties protected by 
such invariants? These two questions are of absolute importance for most of the 
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newly envisioned technological applications of these materials. In the same time, these 
questions are extremely difficult to deal with both analytically and numerically. Progress 
in this direction, even if minuscule, is extremely valuable at this stage. 

We showcased two applications of the Non-Commutative Theory of the Chern 
invariant, which gave us a glimpse of the methods of Non-Commutative Geometry 
and what can be achieved by such methods. For Chern insulators, the analysis gave 
quite explicit and satisfactory answers: we could spell out the most general conditions 
that assures the quantization and invariance of the non-commutative Chern number, 
we were able to conclude the existence of protected extended bulk states and even to 
draw the phase diagram of our simple Chern insulator model. For Quantum spin-Hall 
insulators, we defined the non-commutative spin-Chern number and again were able 
to spell out the most general conditions that assures its quantization and invariance, 
even if in a formal way, and to give some simple explicit physical situations, the weak 
disorder regime, where the quantization and invariance of the spin-Chern number can be 
firmly established. We have argued that the non-commutative spin-Chern number is an 
effective tool for mapping the topological phase of the model with strong disorder. The 
numerical simulations indicated that the spin-Chen number also protects the extended 
bulk states of the QSH insulator. 

The construction of the non-commutative spin-Chern number can be repeated for 
many other situations. There is a firm belief among the condensed matter scientists 
that every discrete symmetry has an associated topological insulator and a topological 
invariant that classifies its topological phases and the efforts towards identifying new 
classes of topological insulators are already underway [1221 11231 1124j . We believe that 
novel topological phases can be found by looking for analogs of a Z) that is, w operators 
such that the spectrum of PwP splits in distinct spectral islands that carry non-zero 
Chern numbers. This is our primary reason to believe that the Non-Commutative 
Theory of the Chern number can become increasingly relevant for the field of topological 
insulators. 
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